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ABSTRACT 

Optimum  receivers  for  detecting  binary  signals  in  additive 
colored  Gaussian  noise  are  analyzed  and  their  performance  evalu- 
ated in  terms  of  bit  error  probabilities  (P  ) .   Implementation 
and  practical  design  implications  of  such  receivers  is  discussed, 

Evaluation  of  P   for  receivers  that  are  optimum  for  additive 

e  r 

white  Gaussian  noise  (WGN)  environments  but  due  to  jamming  or 
"friendly"  ECM  interf erers ,  must  operate  in  a  colored  Gaussian 
noise  environment  has  been  carried  out.   It  was  generally  found 
that  such  receivers  do  not  perform  significantly  worse  than 
receivers  specifically  designed  to  operate  in  a  colored  noise 
environment.   Examples  were  considered  in  which  the  colored 
noise  interference  was  modeled  as  the  output  of  a  one-pole 
filter  driven  by  WGN.   Additional  work  has  been  carried  out  on 
the  jamming  of  binary  (colored  noise)  receivers  using  a 
deterministic  jammer  model.   While  this  modeling  assumption 
needs  to  be  refined,  it  has  been  demonstrated  that  a  power  con- 
strained jammer  can  seriously  degrade  the  performance  of  a 
receiver  designed  to  operate  in  a  colored  noise  environment. 
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I.   INTRODUCTION 

The  theory  of  statistical  signal  detection  and  estimation 
in  the  presence  of  colored  noise  is  described  in  many  text- 
books [Refs.  1,2,  3].   However,  applications  and  practical 
system  design  considerations  as  well  as  implementations  based 
on  the  developed  theory  are  not  often  discussed. 

There  are  some  [Refs.  4,5,6]  DOD  research  publications 
which  deal  with  signal  reception  in  colored  noise.   None  of 
these  publications  uses  colored  noise  thoery  as  developed  in 
textbooks  or  analyzes  advantages  and  disadvantages  of  signal 
designs  that  account  for  the  presence  of  colored  noise. 

The  goals  of  this  thesis  are  to: 

1.  Discuss  some  practical  applications  that  can  be 
derived  from  colored  noise  signal  detection  theory. 

2.  Analyze  practical  design  implications  of  the  theory. 

3.  Present  advantages  as  well  as  disadvantages  of  using 
theoretical  results  involving  colored  noise  inter- 
ferences as  compared  to  results  dealing  with  white 
noise  interference  models. 

Specifically,  the  following  problems  will  be  analyzed: 

1.  The  design  of  a  binary  communication  receiver  in  the 
presence  of  colored  noise  interference  when  the  signals 
used  to  transmit  the  binary  information  are  completely 
known . 

2.  Design  of  an  optimum  signal  set  for  a  receiver  opera- 
tion under  the  same  conditions  as  in  1.  above. 

3.  Comparison  of  the  performance  of  the  receivers  analyzed 
in  1.  and  2.  above  relative  to  an  equivalent  receiver 
designed  to  operate  in  white  noise  interference. 
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4.  Evaluation  of  the  effect  of  an  RF  prefilter  on  the 
performance  of  a  binary  communication  receiver. 

5.  Evaluation  of  the  performance  of  receivers  designed 
for  colored  noise  interference  operating  in  the 
presence  of  jammers. 

This  thesis  is  divided  up  as  follows. 

In  Chapter  II  we  present  briefly  colored  noise  theory  and 
the  integral  equations  governing  the  receiver  design.   In 
Chapter  III  we  analyze  Fredholm  Integral  Equations  and  the 
techniques  used  for  solving  them  for  baseband  signals  and 
bandpass  signals.   In  Chapter  IV  we  discuss  receiver  design 
and  performance  in  the  presence  of  colored  noise  interference, 
according  to  the  rules  of  colored  noise  theory.   We  also  com- 
pare these  results  to  the  more  conventional  receiver  designed 
to  operate  in  a  white  noise  only  environment.   In  Chapter  V 
we  analyze  the  effect  of  using  RF-preamplif iers  in  digital 
receivers.   In  Chapter  VI  the  sensitivity  of  this  receiver  is 
evaluated  when  operating  in  the  presence  of  colored  noise 
interference  and  a  deterministic  jammer  signal  which  is  opti- 
mum in  a  specific  sense.   Performance  comparisons  to  equivalent 
white  noise  environment  receivers  are  presented.   The  Conclusions 
and  interpretations  of  the  results  obtained  are  presented  in 
Chapter  VII. 
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II.   THEORY  OF  COHERENT  SIGNAL  DETECTION  IN  THE 
PRESENCE  OF  COLORED  NOISE 

A.   DETECTION  IN  THE  PRESENCE  OF  WHITE  GAUSSIAN  NOISE 

The  design  of  coherent  signal  receivers  in  the  presence 
of  additive  white  Gaussian  noise  (WGN)  is  widely  dealt  with 
in  the  literature  [Refs.  7,8,9,10] .   The  optimum  receiver 
(in  the  sense  of  producing  minimum  probability  of  error)  for 
discriminating  between  two  different  yet  completely  known 
signals  in  additive  WGN  is  shown  in  Fig.  2.1. 

The  receiver  of  Fig.  2.1  is  optimum  (i.e.,  minimum  proba- 
bility of  error)  when  the  received  signal  z (t)  is  either 

Hypothesis  1^:    z(t)   =   y-,(t)  +  n(t)     0  £  t  £  T 
Hypothesis  HQ :    z(t)   =   yQ  (t)  +  n(t)     0  <_  t   <_   T 

where  y,  and  y_ (t)  are  known  deterministic  signals  and  n(t) 
is  a  sample  function  of  a  WGN  process.   For  convenience  we 
define 


Yd(t)   =   yx(t)  -  y0(t) 


The  threshold  level  shown  in  Fig.  2.1  is  given  by  [Pefs 
7,8,9,10] . 


N 
Y   =  -j   in   I±Z£1.  +  |(E1  -EQ)  (2.1) 
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where 


T 
E.   =    /   y2(t)dt      i  =  0,1 
1      0     X 

N0 
and  -=-  is  the  two-sided  power  spectral  density  level  of  the 

WGN  interference.   Also,  P  is  the  prior  probability  that  signal 

y, (t)  was  sent. 

If  equiprobably  and  equal  energy  signals  are  transmitted, 
then  P  =  1/2  and  the  threshold  y   becomes  zero. 

Assuming  equiprobable  signals,  the  performance  of  this 
receiver  is  given  by 


/E(l-p) 

e  -«--*  x  y        N 


P„       =      ERFCi({\rKt         )  (2.1A) 


0 


where 


T 
=  \     j     tY^t)  +  y^(t)]dt   =   AVG  Bit  Energy 


and 


1    T 
p      =      -      j      y0(t)Yl(t)dt 


0 

The  complementary  function  ERFC* ( . )  used  throughout  this  thesis 
is  defined  as 


00        2  . 
ERFC*(v)   =    /   —  e"X  /2  dx 
v    /2rf 
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For  the  so-called  antipodal  signals,  y1(t)  and  yQ  (t)  are  related 
by 


Yl(t)  =  -y0(t; 


Hence  p  =  -1  and  Eq.  2 . 1A  becomes 


Pe   =   ERFC*(/(2E)/NQ)  (2.2) 


It  is  important  to  notice  that  P   is  independent  of  the 
particular  waveform  shapes  used.   Equation  (2.1A)  demonstrates 
that  the  signal-to-noise  ratio  (SNR)  E/N   and  the  normalized  sig- 
nal correlation  coefficient  p  are  the  only  factors  affecting  P  . 
Such  will  not  be  the  case  when  the  noise  interference  is  colored 


B.   DETECTION  IN  THE  PRESENCE  OF  COLORED  NOISE 

In  certain  cases,  the  transmitted  signals  can  encounter  a 
nonwhite  colored  Gaussian  interference.   The  most  common  such 
cases  arise  when: 

1.  Between  the  actual  white  noise  source  and  the  signal 
processing  part  of  the  receiver,  there  are  some 
bandpass  elements  such  as  antennas  or  RF  filters 
which  shape  the  noise  spectrum  so  that  it  no  longer 
is  white. 

2.  In  addition  to  the  desired  signal  at  the  front  end  of 
the  receiver,  there  is  an  interfering  signal  which 

may  be  some  ECM  jammer  or  may  be  a  "friendly"  electronic 
emitter  causing  interference  in  the  communication 
channel.   In  radar/sonar  systems  such  interference  is 
frequently  caused  by  multiple  targets. 

3.  Multipath  channel  interferences  arise  which  effectively 
add  a  colored  noise  component  to  the  channel. 
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The  basic  decision  model  can  now  be  specified  as  follows 

Hypothesis  H,  :   z  (t)   =   y,  (t)  +  n  (t)  +  n  (t)     0  <  t  <  T 
J-  _I_  w  c  ~~~    — 


Hypothesis  H  :   z  (t)   =   yQ  (t )  +  n  (t)  +  n  (t)     0  <_   t  <_   T 


where  n  (t)  is  the  colored  noise  component  and  n  (t)  is  the 

V—-  vv 

white  noise  component.   Notice  that  a  white  noise  component  is 
present  in  the  model.   It  is  appropriate  to  assume  that  the 
interference  contains  also  an  independent  white  component  due 
to  the  fact  that: 

1.  Practical  systems  always  will  contain  a  nonzero 
thermal  white  noise  component.   Even  shot  noise  which 
is  dominant  in  the  optical  range  of  the  spectrum  is 
also  practically  a  white  noise. 

2.  As  will  be  discussed  in  Appendix  A,  the  white  noise 
component  enables  us  to  guarantee  that  our  mathematical 
solutions  will  be  meaningful. 

The  conventional  approach  in  the  design  of  an  optimum  receiver 
is  to  take  "samples"  of  the  received  signal,  express  the  joint 
probability  density  function  of  these  samples  and  then  to 
determine  the  limiting  form  as  the  samples  are  taken  closer 
together  and  their  number  increases  to  infinity.  These  opera- 
tions become  more  difficult  in  the  case  of  colored  noise  since 
the  samples  may  no  longer  be  statistically  independent  [Refs. 

1,2,10]  • 

The  approaches  taken  when  colored  noise  interference  is 

present  are  to: 

1.   Introduce  a  "Whitening"  filter  to  transform  the  inter- 
ference into  a  white  Gaussian  noise  so  that  use  of 
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the  white  Gaussian  interference  analysis  to  solve  the 
problem  is  possible  [Refs.  1,2]. 

2.   Use  the  Karhunen-Loeve  expansion  [Ref.  3]  . 

The  advantage  of  the  Karhunen-Loeve  expansion  is  that  it 
leads  to  a  series  of  elements,  the  coefficients  of  which  are 
uncorrelated.   These  coefficients  represent  the  signal  "samples" 
or  components  along  the  dimensions  associated  with  each  eigen- 
f unction  g . (t) . 

Clearly  both  the  "Whitening  Filter"  approach  and  the  Karhunen- 
Loeve  expansion  approach  lead  to  exactly  the  same  results. 

C.   DERIVATION  OF  THE  OPTIMUM  RECEIVER  IN  COLORED  NOISE  VIA 
THE  KARHUNEN-LOEVE  EXPANSION  METHOD 

Let  us  assume  that  the  noise  is  colored  with  covariance 

function  K  (t,u).   We  expand  the  noise  covariance  function  in  terms 

of  a  set  of  orthogonal  functions.   We  use  the  Karhunen-Loeve 

expansion  in  which  the  orthogonal  functions  are  the  eigen- 

functions  of  the  integral  equation 


T 
/   Kv(t,u)gi(u)du   =   A^Ct)    i  =  1,2,...  (2.3) 


We  are  now  able  to  expand  the  received  signal  z(t),  the  signals 
y.(t),  i  =  0,1,  and  the  noise  n(t)  in  the  coordinate  system 
specified  by  the  set  {g.(t)}.   That  is 


K 
z(t)   =   fc.i.m  I      z.g. (t)  (2.4) 

K+~   i=l   x  x 
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Figure  2.2    Coherent  Colored  Noise  Correlator  Receiver 
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K 
n(t)   =   £.i.m    I      n.g. (t)  (2.5) 

K+°o    i=l   1  x 


K 

I 


y.(t)   =   £.i.m  I      y.-g.(t)     j  =  0,l  (2.6) 


where 


T 
z.   =    /   z(t)g.(t)dt     i  =  1,2,...  (2.7) 

1      0         1 


T 
ni   =    /   n(t)gi(t)dt     i  =  1,2,...  (2.8) 


T 
y..   =    /   y.(t)g.(t)dt    1=1,2,...  (2.9) 

:1      0     :     X         j  =  0,1 


It  is  reasonable  to  assume  that  the  noise  is  zero  mean.   Then 


E(n. )   =   0  i  =  1,2,  ...  (2.10) 


The  covariance  of  n.  and  n.  will  become 

i       3 


T  T 

Efn.n.)   =   E[  /   n(t)g.(t)dt-  /   n(u)g.(u)du] 
1  :         0         1       0         : 


T  T 

/   g.(t)   /   K  (t,u)g.(u)  du  dt  (2.11) 
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Substituting  Eq.  (2.3)  into  (2.11)  yields 

E(n±n  )   =   X±6  (2.12) 

where  6. .  is  the  Kroenecker  delta.   Equation  (2.12)  demonstrates 
that  the  Karhunen-Loeve  expansion  leads  to  noise  "samples" 
that  are  statistically  independent  since  the  noise  has  been 
assumed  Gaussian. 

We  can  now  express  the  likelihood  ratio  test  involving  the 
K  signal  "samples"  as  follows  [Refs.  4,5,6] 

n     1       i      1  (Zi"yli)2,   n 
II   — — — -  exp{-  T  ^r }   H 

i=l  /2¥TT      l         l Ai      / 

A(zk(t))   =   1 2~  \        Y  (2.13) 

*     i        ,1  (2i-y0i)  ,   Hn 

n   — — — —  expl-  T  x^ }    0 

i=l  /2?X7      ^    2Ai 

i 

where  y    is  the  threshold  defined  by  Eq.  (2.1).   Cancelling 
common  terms,  taking  the  logarithm  and  letting  K  to  infinity 
yields  the  decision  rule 


oo    Z 

.1  t . . .  „. 

1=1   1  1=1   1 


oo     2  °° 

iB(AU(t)»    =     i   ji(yu-yoi)  - 1  I   rr'yoi-yii' 


Hi 

J    Iny  (2.14) 

H0 


If  equiprobable  antipodal  signals  are  used,  then  in   y   =  0 . 
Expressing  the  likelihood  ratio  in  terms  of  the  signals  z(t), 
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yx(t),  and  yQ(t)  by  using  Equations  (2.8),  (2.9),  (2.10), 
respectively,  yields 

T    T  °o   g.  (t)g.  (u) 

/    /  [y,  (u)-yn(u)]z(t)  I      -± 5-1 du  dt 

0    0     X      U  i=l      Ai 

i    T    T  -   g  (t)g. (u)        X 

+  2   /   /  [y0(t)y0(u)-yi(t)yi(u)]  I  r  dudt  >   0 

z    0    0     U     U      X     L  i=l      Ai  < 

Ho 

(2.15) 


We  define 


A  T  g. (t)g. (u) 

h1(t)  -h  (t)   =  h  (t)   =   /  (y1(u)-y()(u))  \     — — r-i du 

0  i=l      i 


(2.16) 


Substituting  Eq.  (2.16)  in  Eq.  (2.15)  yields  the  simplified 
decision  rule 

Hl 
T  ,    T 

/   z(t)h  (t)dt^   i   /  [y  (t)h  (t)-y  (t)h  (t)]dt   =   y'    (2.17) 
0         Q      „     0  uu 

H0 

where  hd(t)  is  defined  byEq.  (2.16).   The  term  on  the  right-hand 
side  of  the  inequality  sign  is  a  constant  and  may  be  considered 
as  a  new  threshold  y ' . 

We  can  get  a  different  mathematical  form  for  h,  (t)  by 
multiplying  Eq.  (2.16)  by  K  (t,u)  and  then  integrating  over 
the  interval  (0,T).   We  thus  obtain 
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T  T    T 

/   Kv(t,u)hd(u)du       /    /   Kv(tfu)  (y]_(u)  -yQ(u)) 


0    0 


00   9^  (t)g.  (u) 

V   _i : 

■  S      A. 
1=1       1 


I        X   x  X       -dtdu 


=  yx(t)  -  yQ(t)  .  (2.18) 


Therefore  h ,  (u)  is  now  defined  by  the  integral  equation 


T 
/   Kv(t,u)hd(u)du   =   y1(t)-yQ(t)  (2.19) 


An  optimum  receiver  structure  can  now  be  obtained  as  a  direct 
consequence  of  the  decision  rule  of  Eq.  (2.17).   Eq.  (2.19) 
defines  h,  (t)  implicitly  and  the  receiver  structure  is  shown 
in  Fig.  2.2.   We  refer  to  this  receiver  as  a  "Colored  Noise 
Receiver."   Fig.  2.2  shows  that  this  receiver  is  a  correlation 
detection  receiver  much  like  a  "White  Noise  Receiver,"  except 
that  the  correlating  signal  is  no  longer  of  the  same  form  as 
that  of  the  transmitted  signal  but  instead  is  given  by  the 
solution  to  the  integral  equation  of  Eq.  (2.19). 

D.   RECEIVER  PERFORMANCE 

The  performance  of  a  binary  communication  receiver  can  be 
quantified  as  the  probability  of  making  an  erroneous  decision 
labeled  P  .   This  involves  finding  the  probabilities  that  the 
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output,  I,  of  the  correlator  (Fig.  2.2)  exceeds  or  is  exceeded 
by  the  threshold  given  knowledge  of  which  signal  was  trans- 
mitted.  This  can  be  accomplished  since  I  is  a  linear  combina- 
tion of  a  Gaussian  noise  and  therefore  is  a  Gaussian  random 
variable.   Thus,  P   is  given  by 


Pe   =   p{H0}P{I>y'|H0}  +  P{H1}P{I<y|H1> 


Y' 


=   P(HQ}    /   P{l|HQ}dI  +  P{H1>    /    P{I|H  }dl 


where  P{Hn>,  P{H1 }  are  the  prior  probabilities  of  sending  signal 

y   or  y, ,  respectively.   Since  we  assume  equiprobable  antipodal 

signals,  P{H_}  =  P{H, }  =  y  and  y'  as  defined  by  Eq.  (2.17)  is 

zero.   Since  I | H_  and  l|H,  are  Gaussian  variables,  their  mean 

and  variance  only  need  to  be  found  in  order  to  evaluate  the 

P  .   The  mean  value  of  I  I Hn    is  given  by 
e  '  0     ^       J 


T  T 

E{I|HQ}   =   E{  /  [yQ(t)+n(t)]hd(t)dt}  =    /   yQ(t)hd(t)dt 


=   mA.  (2.21) 


0 


Similarly 


T 
E{l|H,}   =    /   y, (t)h,(t)dt   =   m. .  (2.21A) 

1      0    1         d         x 
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For  antipodal  signals,  that  is  y. (t)  =  -y, (t) ,  we  have 


ml   =   ~  m0  '  (2.22) 


It  can  be  easily  shown  that  the  variance  of  random  variable 
I  conditioned  by  Hn  or  H,  is  given  by 


T 

Var{l|H  }   =   Var{l|H  }   =   E{ [  /   n (t) h  (t )dt] 2 }  =  a2    (2.23) 
x  u  0         a  I 


Eq.  (2.23)  can  be  written  in  the  form 


T    T 

a2.   =    /    /   h,  (t)K  (t,u)h,(u)  dudt  (2.23A) 

I      Q>       QJ         d     v       d 


Substituting  Eq.  (2.19)  in  Eq .  (2.23A)  yields 


T 
a2   =    /   hd(t)yd(t)  dt   =   21^  (2.24) 


From  Eq.  (2.21)  and  Eq.  (2.24)  and  since  y '  =  0  for  antipodal 
equiprobable  signals,  P   becomes 


-,         °°  -,  (I-m,) 

Pq      =      i      /  x  exp{ ± — }  dl 

0  x/o      2  2aT 

v27raT  I 


I 

1  °  1  d-mn)2 

+  i        /  x  exp{ 1 — }  dl  (2.25 

-00  x/o      2  2aT 

v27Ta  I 
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A  change  of  variables  in  Eq.  (2.25)  and  use  of  Eqs.  (2.21-2.23) 
yields 

Pe   =   ERFC*(  yj  /  yd(t)hd(t)dt)   =   ERFC*(^|.2m1)    (2.26) 

Observe  that  the  performance,  of  the  "Colored  noise"  receiver 
does  depend  on  the  signal  waveforms.   For  white  noise  inter- 
ference Pe  was  shown  to  be  independent  of  the  signal  waveforms. 
(See  Eq.  (2.2).) 

E.   OPTIMUM  SIGNAL  DESIGN 

Since  the  performance  of  the  receiver  analyzed  in  the 
previous  sections  depends  on  y,(t),  there  may  be  an  optimum 
waveform  set  for  minimum  probability  of  error.   From  Eq.  (2.26) 
it  is  clear  that  by  making  y,  (t)  large,  P   can  be  made  small. 
Thus,  to  make  the  optimization  problem  meaningful  an  energy 
constraint  is  placed  on  the  signal  set.   That  is,  with  fixed 


T 

E   =   i  /  (yn(t)  +  y*(t))dt 
c    Q     u 


J   = 


T  T 

/   y  (t)h  (t)dt  -  y[  /  (yj;(t)+y^(t))dt  -  2E)         (2.27) 
0     a    a  0     u      -l 


is  to  be  maximized,  where  u  is  the  Lagrange  multiplier.   It 
then  follows  from  the  calculus  of  variations  [Refs.  1,2]  that 
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the  optimum  signal  set  obeys 

y1(t)   =   -  yQ(t)  (2.28) 


and 


T 
/   Kv(t,u)y1(u)du   =   Aiyi(t)  (2.29) 


or  equivalently 


T 
/   Kv(t/U)yd(u)du   =   A±yd(t)  (2.30) 


There  are  many  solutions  to  Eq.  (2.30)  and  the  one  which 
corresponds  to  the  minimal  A.  should  be  chosen  [Refs.  2,14]. 

F.   SUMMARY 

This  chapter  discusses  the  theory  of  designing  a  receiver 
in  the  presence  of  colored  noise  interference.   Figure  2.2 
shows  the  block  diagram  of  this  receiver  and  reveals  the  fact 
that  this  receiver  is  basically  a  correlation  detection  re- 
ceiver.  The  only  difference  is  that  the  correlating  signal 
h, (t)  which  is  given  by  the  solution  of  the  integral  equation 


0 


T 
/   Kv(t,u)hd(u)du   =   Yl(t)  -  yQ(t)   =   yd(t)    (2.19) 
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must  be  used  in  place  of  y,  (t)  as  is  done  for  an  optimum 
receiver  operating  in  WGN  interference.   The  performance  of 
this  receiver  assuming  equiprobable  antipodal  signals  is  given 
by  equation 


T 


=   ERFC*(\/j   /   y,(t)h,(t)dt)  (2.26) 


0 


The  performance  depends  heavily  on  the  signal  waveforms. 
An  optimum  signal  waveform  set  is  given  by  the  solution  of 
the  integral  equation 


/   Kv(t,u)yd(u)du   =  y~   yd(t)  .  (2.30) 

0  i 


The  design  procedure  of  a  colored  noise  receiver  will  then 
consist  of  the  following  steps: 

1.  Identifying  the  environmental  noise  and  formulating 
its  correlation  function. 

2.  Solving  Eq.  (2.19)  for  the  correlating  signal  h ,  (t) 
when  y. (t)  and  yn (t)  are  known. 

3.  If  y^(t)  or  yg(t)  are  not  given,  then  Eq.  (2.30)  must 
be  solved  first  for  an  optimal  signal  set  and  only 
then  can  Eq.  (2.19)  be  solved  for  the  optimum 
correlating  signal. 
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III.   SOLUTION  OF  FREDHOLM  INTEGRAL  EQUATIONS 

A.   FREDHOLM  INTEGRAL  EQUATIONS 

As  demonstrated  in  the  previous  chapter,  an  important 
step  in  the  design  of  the  colored  noise  receiver  involves 
solving  an  integral  equation  in  order  to  obtain  h ,  (t) .   That 
is,  a  solution  to  the  following  equation  must  be  found. 


T 

/   K  (t,u)h,(u)du   =   y,(t)     0  <  t  <  T  (2.19) 

'        v      d  -^d         —   — 


This  equation  is  called  a  Fredholm  Equation  of  the  First 
Kind.   The  function  K  (t,u),  namely  the  noise  covariance,  is 
called  the  kernel  of  the  equation.   If  the  kernel  of  Eq.  (2.19) 
contains  singularities,  or  equivalently  if  the  colored  noise 
contains  an  additive  white  noise  component,  then  K  (t,u) 
takes  on  the  mathematical  form 

Nl 
Kv(t,u)   =   -y6(t-u)  +  Kc(t,u)  (*)(3.1) 

which  when  substituted  in  Eq.  (2.19),  yields 


N,  T 

-^h.(t)  +   /   K  (t,u)h,(u)du  =   y,  (t)   0  <  t    T     (3.2) 
Z      a.  ~J         c  d  d       —   — 


( *) 

v  'Throughout  this  thesis,  6(t)  denotes  the  Dirac  Delta 

Function. 


29 


This  equation  is  called  a  Fredholm  Equation  of  the  Second 
Kind. 

The  properties  satisfied  by  these  equations  have  been  dis- 
cussed and  proved  in  many  textbooks  [Refs.  2,11,12].   We  shall 
state  here  only  those  properties  that  are  important  to  the 
present  work. 

Property  1:   If  the  kernel  does  not  contain  singularities 
(i.e.,  no  white  noise  component)  a  finite  square  integrable 
solution  to  the  Fredholm  I  equation  will  not  exist. 

Property  2:   In  this  case  of  kernel  singularities,  a 
solution  to  the  Fredholm  I  equation  will  exist  only  if  we 
allow  it  to  contain  singularity  functions  (impulses). 

The  solution  will  then  be  of  the  form 


h,(t)   =   h  (t)  +  I    a.h,   (t)  +  I   b  6(k)(t)  (3.3) 

d         p      V   l  hi      k   k 


where  h  (t)  and  h,  . (t)  are  the  particular  and  homogeneous 

solutions,  respectively,  to  a  differential  equation  derived 

(k) 
from  the  Fredholm  I  equation  and  6    (t)  is  the  k-th  derivative 

of  6  (t)  [Ref .  1] . 

Property  3:   The  solution  to  Fredholm  equations  are  at 

best  tedious  to  obtain  and  in  many  cases  solutions  are  very 

difficult  or  impossible  to  obtain.   In  two  specific  cases 

there  is  a  straightforward  procedure  for  solving  Fredholm  II 

equations . 
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1.  If  the  kernel  is  separable,  the  solution  is  quite 
easy  to  obtain  [Ref.  11]. 

2.  If  the  noise  power  spectral  density  is  a  ratio 
of  two  polynomials,  a  solution  can  be  obtained 
after  following  a  specific  procedure.   This  situa- 
tion occurs  when  the  colored  noise  is  the  steady 
state  response  of  a  linear  time  invariant  system 
excited  by  white  noise. 

In  this  research,  we  shall  deal  only  with  Fredholm  II  equations 
for  the  following  reasons. 

1.  From  a  practical  standpoint,  we  do  not  want  to  deal 
with  the  problem  of  trying  to  generate  impulse  functions. 

2.  In  real  physical  systems  there  will  always  be  some 
white  noise  component,  however  small,  due  to  thermal 
effects  in  the  electronic  circuitry.   One  is  never  able 
to  totally  eliminate  the  white  noise  component. 

Also,  we  will  deal  only  with  colored  noise  having  a  rational 

spectra  since  it  best  models  the  output  of  real  physical 

systems. 


B.   GENERAL  SOLUTION  TO  FREDHOLM  II  EQUATIONS  FOR  BASEBAND 
SIGNALS 


The  Fredholm  II  equation  of  interest  is 


N,  T 

^h,(t)  +   /   K  (t,u)h  (u)du   =   y  (t)   0  <  t  <  T     (3.2) 


We  assume  that: 

1.   The  white  noise  component  has  power  spectral  density 
level  N./2  watts/Hz. 
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2.  The  noises  are  wide-sense  stationary  (W.S.S.) 

3.  The  Power  Spectral  Density  (P.S.D.)  of  the  colored  noise 
can  be  expressed  as  a  ratio  of  two  polynomials  that  are 
functions  of  the  complex  variables,  namely 


,  ,  .      N(s2) 
D(sZ) 


(<f)c(s)  and  K-Ct)  are  a  two  sided  Laplace  transform 
pair) .   Multiplying  both  sides  of  the  above  equation  by 
D(s2)  yields 

D(s2)cJ>  (s)   =   N(s2)  (3.4) 


Multiplication  by  s  corresponds  to  differentiation  with  respect 
to  t  in  the  time  domain.   So,  Eq.  (3.4)  becomes 


D(p2)K  (t-u)   =   N(p2)6(t-u)  (3.5) 

c 


2 

where  p  =  d/dt.   Operating  on  Eq.  (3.2)  with  D(p  )  yields 


N  T 

D(p2)  [-i  hd(t)]  +   /   hd(u)D(p2)  [Kc(t,u)]du 


=   D(p2)  [y,(t)]  (3.6) 


Substituting  Eq.  (3.5)  in  (3.6)  and  performing  the  integration 
yields 


D(p2)  [-j  hd(t)]  +  N(p2)  [hd(t)]   =   D(p2)yd(t)  (3.7) 
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Eq.  (3.7)  is  a  differential  equation  that  must  be  solved 
completely.   In  this  research  we  deal  with  a  simple  colored 
noise  model  in  which  the  output  of  a  1st  order  Butterworth 
filter  driven  by  white  noise  is  taken  as  the   source  of  colored 
noise.   Thus  its  P.S.D.  is  given  by 


2 

cj>  (s)   =   5—  =   — ^ j  (3.8) 

D(s^)      -sz  +  3 


Equivalently  in  the  time  domain 


Kc(x)   =   ae  3 I T I  (3.9) 


Substituting   Eq.     (3.8)    into    (3.7)    yields 


N 
D(p2)  [^  hd(t)]    +   2a3[hd(t)]      =      D[p2]yd(t) 


2 
Operating  with   D(p    )    yields 


N      ..  N  "2 

-±  h    (t)    +    (^3      +   2aB)hd(t)      =     -yd(t)   +3   yd(t)    .       (3.10) 


Eq.  (3.10)  is  a  second  order  differential  equation.   Its 
solution  is  of  the  form 


h(t)   =  h  (t)  +  K1hhl(t)  +  K2hh2(t)  (3.11) 


where  h  (t)  is  the  particular  solution  and  h,  ,(t)  an^  h,  2 (t) 
form  the  homogeneous  solution.   Substitution  of  Eq.  (3.11) 
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into  Eq.  (3.2)  leads  to  two  simultaneous  equations  that  K, 
and  K~  must  satisfy.   Solving  for  K,  and  K^  explicitly  gives 
the  complete  solution. 

A  similar  procedure  must  be  applied  when  the  noise  is 
modeled  as  the  output  of  a  higher  order  Butterworth  filter. 
For  such  a  case  (Nth  order  Butterworth  filter) ,  the  noise 
P.S.D.  is 


Q2N  XT/  2^ 

6  (s)   =   2 =   N(s  }  (3  12) 

(3)    +s       D(s  ) 


Substituting  Eq.  (3.12)  in  Eq.  (3.2)  yields 


N    2N  N 

-^[^-2-hd(t)]  +  [^(e)2N  +  (6)2N]hd(t) 

.2N 
=   ^2N^d(t)  f  ^)2yd(t)  (3.13) 


There  will  now  be  2N  homogeneous  solutions  and  a  particular 
solution.   Substitution  into  the  Fredholm  II  equation  will 
lead  to  2N  simultaneous  equations  from  which  the  constants 
associated  with  the  homogeneous  solutions  must  be  determined. 

The  complete  solution  for  the  case  in  which  y,  (t)  is 
rectangular  or  sinusoidal  is  worked  out  in  detail  in  Appendices 
A  and  B.   The  procedure  is  long  and  tedious  so  that  for  higher 
order  filters  numerical  techniques  must  be  utilized. 
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C.   GENERAL  SOLUTION  TO  FREDHOLM  II  EQUATIONS  FOR 
BANDPASS  SIGNALS 

In  practice, the  communication  signals  are  baseband  signals 

which  modulate  a  carrier  prior  to  transmission.   The  received 

signals  can  be  modeled  as 

Hypothesis  H,  :   z(t)   =   y,  (t)cos  u)Qt  +  n(t)   0  <_   t  <_  T 

(3.14) 
Hypothesis  HQ :   z(t)   =   yQ(t)cos  wQt  +  n(t)   0  <  t  <_   T 

yd(t)  cos  wQt   =   [y1(t)-yQ  (t)  ]cos  a)Qt  (3.15) 

where  oo^  is  the  carrier  frequency  and  y,  (t)  ,  and  yn(t)  are 
the  baseband  signals. 

The  Fredholm  II  equation  now  becomes 


N,  ~  T  ~ 

-^h,(t)  +   /   K  (t,u)h,(u)du   =   y-,(t)  cosoi.t         (3.16) 
2d      ())         cd  J  d  0 


where  h(t) ,  K(t-u)  represent  bandpass  waveforms. 

The  bandpass  autocorrelation  function  K  (t,u)  can  be  ex- 
pressed as 


K  (t-u)   =   K  (t-u)cos  a)n(t-u)  (3.17) 

C  CO 


where  K  (t-u)  is  the  baseband  autocorrelation  function, 
c 
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The  solution  of  this  equation  follows  the  same  procedure 
as  the  one  used  in  the  solution  of  Eq.  (3.2)  but  is  consider- 
ably more  tedious. 

In  Appendix  C,  we  prove  that  if  the  carrier  frequency  con 
is  much  bigger  than  the  bandwidth  of  the  noise  or  the  band- 
width of  the  information  signals,  then  the  solution  to  Eq. 
(3.16)  is  approximated  by 


hd(t)   =   hd(t)cos  co0t  (3.18) 


where  h,(t)  is  the  solution  to  the  Fredholm  II  (Baseband) 
a 

equation  of  Eq.  (3.2). 

Since  in  practical  cases  the  carrier  frequency  is  much 
bigger  than  the  bandwidth  of  the  data  or  of  the  colored  noise, 
we  will  only  solve  the  Fredholm  II  equation  for  baseband  signals 
and  will  use  this  solution  as  the  solution  for  the  bandpass 
case  using  Eq.  (3.18). 
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IV.   RECEIVER  DESIGN  AND  PERFORMANCE  IN  THE 
PRESENCE  OF  COLORED  NOISE  INTERFERENCE 

A.   INTRODUCTION 

This  chapter  presents  the  analysis  of  a  typical  case  in 
which  the  communication  receiver  operates  in  an  environment 
which  consists  of  colored  noise  interference.   This  interference 
can  be  due  to  an  ECM  jammer,  a  'friendly1  electronic  emitter, 
or  some  multipath  interference.   This  chapter  presents  the 
design  of  the  colored  noise  receiver  which  is  optimized  to 
the  presence  of  the  interference.   The  performance  of  this 
receiver  is  analyzed  and  compared  to  the  performance  of  a 
coherent  digital  communication  receiver  (designed  for  white 
noise  only  interference)  operating  in  the  same  environment. 
The  most  important  parameter  in  the  design  procedure  is  the 
interference  P.S.D.   This  fact  creates  problems  when  designing 
a  colored  noise  receiver  whose  function  is  to  suppress  hostile 
interferences.   However  there  are  many  applications  in  which 
the  P.S.D.  of  the  colored  noise  is  either  known  or  can  be 
measured  sufficiently  accurately.   A  typical  such  application 
occurs  when  a  digital  communication  receiver  has  to  operate  on 
board  a  ship  or  an  aircraft  in  the  presence  of  other  friendly 
emitters  such  as  radars,  ECM  transmitters,  or  navigational 
equipment.   Those  emitters  whose  characteristics  are  known, 
often  cause  significant  degradation  in  the  quality  of  the 
digital  communication  channel.   Even  multipath  interference  can 
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be  measured  and  modeled  as  colored  noise  interference  whose 
power  spectral  density  or  autocorrelation  function  are  known. 
The  results  of  this  chapter  demonstrate  that  in  such  situations, 
utilization  of  a  colored  noise  receiver  with  "proper"  signal 
waveforms  can  improve  the  P   performance. 

B.   THE  MODEL 

The  system  model  consists  of  a  digital  coherent  communica- 
tion receiver  operating  in  the  presence  of  both  colored  noise 
and  additive  white  noise  interference  as  shown  in  Fig.  4.1. 
Both  noises  are  assumed  Gaussian. 

The  signals  are  binary,  at  baseband,  and  encounter  baseband 
interference.   Extensions  to  bandpass  signal  analysis  is 
straightforward,  given  the  results  described  in  Chapter  III, 
Section  C. 

The  colored  noise  source  block  diagram  is  shown  in  Fig.  4.2. 
This  is  the  typical  block  diagram  of  a  noise  jammer  [Ref .  13] . 

The  output  P.S.D.  of  the  colored  noise  source  is 


N0 
Nc(f)   =  -±   H(f)H*(f)  (4.1) 


where  H(f)  is  the  transfer  function  of  the  amplifier  chain. 
For  the  sake  of  simplicity,  we  model  the  amplifier  chain  as  a 
one-stage  amplifier  of  gain  G  followed  by  a  first  order  Butter- 
worth  filter  with  3  db  bandwidth  26.   If  the  input  to  that 
source  is  a  white  noise  with  P.S.D.  level  N ' /2 ,  then  the  output 
P.S.D.  is 


38 


Colored  Noise 


nc(t) 


Signal 


Y1  or  yQ 


White  Noise 
2 


Figure   4.1        Colored  Noise    Receiver   and   Colored 
Noise    Interference 
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Figure  4.3   White  Noise  Receiver  and  Colored 
Noise  Interference 
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N  *          2  N  2 

,  (S)   =   -°  G2  ^  =   -£     g                      M  2) 

i  (S)                 2  G         2  2    2    ,  (4-2> 

-s  +  3  -s  +  3 


where 


N       N ' 
0       0   2 

-j-  =  "T  G  (4'3) 


The  P.S.D.  can  also  be  written  in  accordance  with  the  nota- 
tion of  Appendix  A  as 


/  n        2a3  ,  A     A\ 

<£c(s)   =   — 22  (4.4) 

-s  +3 


where 


N0 
=  -j-  6.  (4.5) 


In  the  time  domain,  the  autocorrelation  of  the  colored  noise 
is  given  by 


Kc(T)   =   a  exp(-  3  |t|)  .  (4.6) 


and  the  power  of  the  colored  noise  is  given  by 


N0 
Pj   =   Kc(0)   =   a   =   -^  3  (4.7) 


One  should  notice  that  since  a  constant  power  source  is 
assumed,  an  increase  in  the  bandwidth  3  must  be  accompanied  by 
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a  drop  in  the  gain  G  so  that  P.  in  Eq.  (4.7)  will  remain 
unchanged  as  S  is  varied.   The  colored  noise  receiver  struc- 
ture is  shown  in  Fig.  4.1. 

The  input  to  the  receiver  is 

Hypothesis  HQ:     z(t)   =   yQ(t)  +  n  (t)  +  n  (t) 

0  <_  t    <_   T 
Hypothesis  H,  :     z  (t)   =   y-,(t)  +  n  (t)  +  n  (t) 


where  y     and  yn    are  defined  as 


/  -A     0  <_   t  <_  T 
Y0(t)   =  (4.8) 

(  0      0  >  t,  t  >  T 


|  A      0  <_   t  <    T 
y1(t)   =  (4.9) 

'  0      0  >  t,  t  >  T 


and  for  convenience  we  define  (as  before) 


yd(t)   =   Yl(t)  -  yQ(t) 


(4.10) 


C.   RECEIVER  DESIGN  FOR  RECTANGULAR  PULSES 

As  discussed  in  Chapter  II,  the  optimum  receiver  is  (as 
can  be  seen  in  Fig.  4.1)  a  correlation  detection  receiver. 
The  correlating  signal,  h ,  (t)  ,  is  the  solution  of  the 
Fredholm  II  Integral  equation 
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N,  T 

~Y   hd(t)  +   /   Kc(t,u)hd(u)du   =   yd(t)  (3.2) 


The  detailed  solution  of  Eq.  (3.2)  for  the  specific  autocorre- 
lation function  and  signals  given  by  Eqs .  (4.6)  and  (4.8),  (4.9) 
respectively,  is  worked  out  in  Appendix  A.   The  solution  is 
given  by  Eq.  (A. 13),  namely 


h^t)   =   C(l  +  K1eyt    +  K  e  Yt)    0  £  t  £  T        (A. 13) 


where  K^r    K2 ,  C,  and  y    are  constants  defined  by  Eqs.  (A. 15), 
(A. 16),  (A. 17)  respectively.   Defining  now 


mi  -  V '  +  ^  <4-11) 


where  Nn  is  given  by  Eq.  (4.3)  and  N-,/2  is  the  white  noise 
P.S.D.  level  and  also  defining 


E   =   3T  (4.12) 


where  T  is  the  length  of  the  integration  time  in  the  receiver, 
it  is  now  possible  to  specify  h, (t)  in  compact  form.   Observe 
that  the  factor  1/T  can  be  interpreted  as  the  bandwidth  of 
the  receiver,  so  that  E  can  be  viewed  as  the  ratio  of  the 
interference  bandwidth  to  the  receiver  bandwidth.   Substituting 
into  the  results  of  Appendix  A  the  above  definitions  as  well 
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as  the  definition  of  a  as  given  by  Eq.  (4.5) ,  we  can  express 
h,(t)  as  a  function  of  A,  NQ,  N, ,  m, ,  and  E  as  follows 


C   = 


4A 


N0+Nl 


(4.13) 


K. 


-Em,    m,  -1   -2Em, 

m, +1    m,-l   -2Em. 
m, -1    m, +1 


(4.14) 


K. 


m,  +1     -Em. 


m, +1     -2Em,  m,-l 
-  e 


mT-T 


m1+l 


(4.15) 


+3m,  t       -3m,  t 

h,(t)   =  C  +  CK,  e      +  CK„e 
a  l  2. 


(4.16) 


Ignoring  the  constant  of  proportionality  C,  it  is  easy  to 
see  that  h,(t)  is  a  function  of  m,  and  E.   Using  Eqs .  (4.7), 
(4.11) ,  (4.12) ,  m,  can  be  written  as 


m. 


N 


0 


1+*I  - 


I  P.4A  T 

1  +  -1 


6N1A2T 


[^  4(JSR)  (SNR)    (4.17) 


where 


JSR  = 


P  . 

_1 

A2  ' 


SNR   = 


2 

AT 


(4.18) 
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The  factor  JSR  is  the  ratio  of  the  jamming  power  to  signal 
power.   The  factor  SNR  is  the  ratio  of  bit  energy  to  white 
noise  power  spectral  density  level.   This  factor  can  also  be 
interpreted  as  the  signal  to  noise  ratio.   The  product 
(JSR) (SNR)  represents  implicitly  the  ratio  of  the  interference 
power  to  the  white  noise  power  at  the  input  to  the  receiver. 
It  is  independent  of  the  signal  power. 

Figure  4.4  presents  a  plot  of  h, (t)  as  a  function  of  time 
normalized  to  T  for  (JSR) (SNR)  =  1,  and  different  values  of 
E.   Figures  4.5  and  4.6  show  similar  plots  for  different 
(JSR) (SNR)  values  with  E  as  a  parameter. 

These  figures  show  that  decreasing  the  interference  power 
or  increasing  E  tend  to  make  the  colored  noise  interference 
less  dominant  in  comparison  to  the  white  noise.   Effectively 
this  makes  the  receiver  behave  very  much  like  a  white  noise 
receiver.   It  is  not  particularly  difficult  to  design  a  system 
whose  output  will  be  h,  (t).   Such  an  implementation  is  suggested 
in  Appendix  D.   Furthermore,  there  are  now  programmable  signal 
generators  in  the  commercial  market.   However  we  should  be 
aware  that  h  ,  (t)  depends  on: 

1.  The  colored  noise  autocorrelation  function,  and  signal 
waveform. 

2.  The  colored  noise  power  relative  to  the  white  noise 
power . 

3.  The  colored  noise  bandwidth  relative  to  receiver 
bandwidth. 

These  factors  must  not  only  be  known  but  must  also  remain  time 

independent,  unless  adaptive  techniques  are  used.   This  would 
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otherwise  seriously  complicate  the  electronic  circuitry.   These 
h^(t)  dependencies  cause  significant  constraints  on  the  application 
and  use  of  colored  noise  theory  to  receiver  design  problems. 

D.   RECEIVER  DESIGN  FOR  OPTIMAL  PULSES 

Colored  noise  theory  results  demonstrate  that  the  per- 
formance of  the  optimal  colored  noise  receiver,  unlike  the 
white  noise  receiver,  depends  on  signal  waveforms  yn(t)  and 
y,(t).   The  optimal  choice  for  yQ  (t)  and  y-.(t)  is  obtained  as 
a  solution  of  the  integral  equation  given  by  Eq.  (2.30) .   Wi 
now  assume  that  the  model  described  in  Section  IV. B  is  valid 
except  that  the  binary  communication  signals  are  no  longer 
rectangular  pulses  but  can  be  chosen  by  the  system  designer. 
In  other  words  the  system  designer  has  one  more  "degree  of 
freedom."   In  order  to  determine  the  optimum  waveforms  to 
be  utilized  by  the  system  designer,  we  solve  the  integral 
equation  of  Eq.  (2.30)  for  the  case  in  which  the  kernel  is 
defined  by  Eq.  (4.6).   The  solution  is  worked  out  in  detail 
in  [Ref.  14]  and  consists  of  a  set  of  cosines  and  sines  of 
frequencies  b.  which  are  the  solutions 

b. 
(tan  b±T  +  -±) (tan  b±T  -  ^-)       =      0,   i  =  1,2,...        (4.19) 

i 

There  is  an  infinite  number  of  solutions  to  Eq.  (4.19) .   Since 
antipodal  signaling  can  be  shown  to  be  optimum,  we  choose 
one  solution  for  y-,(t)  and  yQ  (t)  given  by 
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Y1{t)       =      A  sinbt 


0  <  t  <  T 


Y0(t)   =   -A  sinbt 


y^(t)   =   2A  sinbt 


0  <  t  <  T 


0  <  t  <  T 


(4.20) 


Yx(t)   =   yQ(t)   =   0 


t  >  T   or   t  <  0 


and  analyze  its  effect  on  receiver  design  and  performance. 
Having  specified  y-j(t),  the  correlating  signal  h,  (t)  must  be 
found  as  a  solution  to 


N, 


h,(t)  +   f   K  (t-u)h,(u)du   =   2A  sinbt 
0  d 


(4.21) 


The  detailed  solution  of  Eq.  (4.21)  is  worked  out  in-  Appendix 
Bv  -  The  solution  is  given  by  Eq.  (B.6)r  namely 


h,(t)   =  O  sinbt  +  CK,  eyt  +  CK0e 
a  1        z 

0  <  t  <  T 


-yt 


(B.6) 


where  K..,  K2,  C,  and  Y1  are  defined  by  Eqs .  (B.8),  (B.9),  (B.5) 
and  (A. 17)  respectively.   With  m,  and  E  defined  by  Eq .  (4.11) 
and  Eq.  (4.12)  respectively,  we  can  express  h, (t)  as  a  function 
of  A,  Nn /  N, ,  m, ,  E,  and  b/3 ,  the  latter  being  the  ratio  of 
the  signal's  frequency  to  the  bandwidth  of  the  interference. 
Thus,  the  constants  of  Eq.  (B.6)  are  given  by 
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4A[(|)2  +1] 

C   = u-j-  (4.22) 

N0+N1+Nl(^) 

,  -Em,   .  -2Em, 

[sin  bT  +  Scos  bT](m.+l)e    -  (§■)  (m,-l)  e 

1  .  2  m,+l    m,-l   -2Em1  \**u) 

[1  +  (£)  I-i-y  -  -i---  e     X] 

p    m,  -1    m,  +1 

h  -Em,  , 

-(§0(m,  +1)  +  (m,-l)e     [sin  bT  +  ^  cos  bT] 

K9   =   — ^Ti oB ; (4.24) 

2  ■  2  ni.+l     2Em,  m,-l 

[1  +   (3-      [- r  -  £ 

p     m, -1         m, +1 


Ignoring  the  constant  of  proportionality  C ,  it  can  be  seen 

that  h,  (t)  is  a  function  of  m,  ,  E,  and  -r-  .   [Notice  that 

a  1  p 

sin  bT  can  be  expressed  as  sin(^-3T)  =  sinC^E).]   The  meaning 

p  p 

of  the  factors  m,  and  E  has  been  discussed  previously.   Fig. 

4.7  is  a  plot  of  h,  (t)  as  a  function  of  time  normalized  to 

T  for  (JSR) (SNR)  =  1,  b/3  =  1  and  various  values  of  E.   Fig. 

4.8  is  a  similar  plot  except  that  (JSR) (SNR)  =  10.   Figs.  4.9- 
4.12  are  repetitions  of  Figs.  4.7  and  4.8  for  various  values 
of  b/3. 


E.   PERFORMANCE  OF  THE  "COLORED  NOISE  RECEIVER" 

Once  the  colored  noise  receiver  has  been  designed,  its 

performance  must  be  evaluated.   This  section  analyzes  the 

performance  of  the  colored  noise  receiver  designed  for  the 

detection  of  rectangular  binary  signals. 

The  performance  of  a  colored  noise  receiver  with  equi- 

probable  binary  antipodal  signals  was  derived  in  Chapter  II. C, 
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and  its  probability  of  error  is  given  by  Eq.  (2.26),  namely 


1   rT  1/2 

=   ERFC*(T   /   y,(t)h,(t)dt) 


(2.26) 


If  we    substitute   Yd(t)    as    defined  by  Eqs .     (4.8)    and    (4.9)    and 
substitute   hd(t)    as   defined  by   Eq.     (4.16),    we   obtain 


-,         T      fl    2  3m   t  -3m   t  1/2 

=      ERFC*(?      /      N°A+N    [1  +Kl£  +K2£  ]dt) 


0  0         1 


=      ERFC. 


( 


2 

2A   T 

N 


V1+n7> 


1    + 


3m   T 
K1(e      L   -1) 

3mxT 


-3m  T 
K2(l  -e         L    ) 

3m1T 


(4.25) 


Substituting   Eqs.     ( 4 . 12) - ( 4 .16)     in   Eq.     (4.25)    yields 


P       =     ERFC 


2(SNR) 


[1+(JSR)  (SNR)i] 


-1  + 


-Em         -2E  (m1  +1) 
2(m1+l)-4e         -2e  (nu-1) 

m,+l       m,-l     -2Em, 

Em,  ( t -zr-  e  ) 

1  m,  -1       m,  +1 


1/2 


(4.26) 


Observe  that  if  E  becomes  unbounded,  or  equivalently  if  the 
colored  noise  has  such  a  large  bandwidth  that  its  P.S.D.  level 
is  nearly  zero  for  all  frequencies,  then  the  colored  noise 
receiver  should  reach  the  performance  of  the  white  noise 
receiver.   Indeed,  letting  E  ■*  °°  in  Eq.  (4.26)  yields 
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1/2 
P    =   ERFC^(2SNR)  '  (4.27) 


where  Eq.  (4.27)  is  the  same  as  Eq.  (2.2)  for  the  white  noise 
receiver. 

Figs.  4.13  and  4.14  show  the  performance  of  this  receiver 

for  SNR  =  10  and  SNR  =  1  respectively,  for  various  values  of 

JSR.   It  can  be  easily  seen  that  as  JSR  increases.  P   increases 

"*  e 

-3 
also.   In  fact,  if  we  define  P   =10    error/bit  as  the  maximum 

e 

probability  of  error  tolerable,  looking  at  Fig.  4.12,  one  can 
say  that  the  receiver  will  not  function  properly  for  JSR 
greater  than  -12  db .   These  graphs  show  also  the  effect  of 
increasing  E.   Increasing  E  spreads  the  jamming  power  over 
larger  frequencies  thus  making  the  P.S.D.  level  lower  at  all 
frequencies.   This  causes  a  decrease  in  the  amount  of  channel 
interference  which  in  turn  causes  an  improvement  in  the 
receiver  performance. 

F.   PERFORMANCE  OF  THE  WHITE  NOISE  RECEIVER 

In  order  to  better  understand  the  performance  of  the 
colored  noise  receiver,  it  is  desirable  to  compare  its  per- 
formance to  that  of  the  coherent  digital  white  noise  receiver 
when  both  operate  under  the  same  conditions.   In  other  words, 
we  deal  with  the  model  described  in  Section  B.   However,  the 
colored  noise  receiver  of  Fig.  4.1  is  replaced  by  the  white 
noise  receiver  of  Fig.  4.3. 
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The  input  to  the  receiver  is 

Hypothesis  H  :     z(t)   =   yQ(t)  +  n  (t)  +  n  (t) 

0  £  t  <_   T 

Hypothesis  H, :     z(t)   =   y, (t)  +  n  (t)  +  n  (t) 
i  -'lew 

where  y~ (t)  and  y, (t)  are  assumed  to  be  two  antipodal  signals. 
The  input  signal  z(t)  is  correlated  with  the  signal  y.(t)/2/E 
where  E   is  the  energy  of  each  signal.   That  is, 


T  T 

E„   =   /   y^(t)dt   =    /   y^(t)dt  (4.28) 

0     l  0     ° 


The  receiver  generates  the  statistic  I ,    which  is  a  Gaussian 
random  variable.   Its  mean  under  both  hypotheses  is  given  by 

1     T 
E{£|H  ]   =   EI-^   /   y,  (t)  [y  (t)+n  (t)+n  (t)]dtj 

±  /E   0 

s 


T 

-±-      j      y^(t)dt   =   /E~  (4.29) 

/E   0     L  S 

s 


E[£|H  ]   =   -  /Es  (4.30) 


Eqs .  (4.29)  and  (4.30)  were  derived  under  the  assumption  that 
both  noises  are  zero  mean.   This  assumption  is  reasonable 
given  the  physical  sources  of  most  noises. 
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The  conditional  variance  of  I    is  given  by 


Var[£|H1]   =   E [ ( l-E ( i | E± ) ) 2 | E±] 


E{  (  /    (y,  (t)+n  (t)+n  (t)dt -/E  at  H   (4.32) 

0     /E     i      c      w  s 

s 


..     T    T 
g*-  /   /  Y1(t)y1(T)  [E{n  (t)n  (t)} 

s  0    0     x     x       c    c 


+  E{n  (t)n  (t)}  +  E{n  (t)n  (x) } 
WW  c     w 


+  E{n  (T)n  (x)  }]dt  dx 
c    w 


Since  n  (t)  and  n  (t)  are  assumed  to  be  statistically  indepen- 
c  w  J  f 

dent  zero  mean  random  processes  ,  we  have 


E[n  (x)n  (t)]   =   E[n  (t)n  (x)]   =   0  (4.33) 

c    w  c    w 


Furthermore 


Nn3      -ft It-   rl 
E[nc(t)nc(x)]       =      -^-EplT:Tl  (4.34) 


and 


N0 
E{n    (x)n    (t)}      =      -^-6(t-x)  (4.35) 

WW  i. 


62 


Substituting   Eqs.     (4.33),     (4.34),     (4.35)     into   Eq .     (4.32) 
yields 

T         T   A2Nnge'e ' t_T'  „    N, 

Var[£|Hn]       =      -L      f         f     [ \i +  AZ   -=£  ( t-x  )  ]  dtdi       (4.36) 

1  s    0         0  q  Z 

Performing   the    integration  yields 


VarU|Hl]      =     ^  +  ^[E-1   Yxp('E)]  (4.37) 


where  E  is  defined  by  Eq.  (4.12). 

Applying  the  same  procedure  for  the  evaluation  of 
Var[£|H_],  we  can  easily  show  that 


VarUlH-jJ   =   VarU|HQ]  .  (4.38) 


Knowing  the  statistical  behavior  of  I,    the  performance  of  the 
receiver  can  now  be  derived.   We  obtain 


P_   =   BBPC^/jj g ^ >  (4.39) 

_i   0  rE-l  +  exp(-E)  , 
2    2  L      E       J 


2E  /N, 

ERFCA  / N   S  L ) 

1    0  fE-l  +exp(-E) 1 
N,  L       E      J 

The  ratio  NQ/N1  can  be  written  (using  Eq.  (4.18))  as 
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Nn      4P.      4P.A2T      .    P.    2        . 

W~      =      6N1   =     J  2    =   E(4)(V")   =   |(JSR)(SNR)       (4.40) 


Substituting  Eqs .  (4.18)  and  (4.40)  into  Eq.  (4.39)  yields 


Pp   =   ERFC*{  / ^SNR }  (4#41) 

VI  +4(JSR)  (SNR)  [E-l  +exp(-E)]/E 

Eq.  (4.41)  specifies  the  performance  of  the  white  noise  receiver 
in  the  presence  of  colored  noise  interference  as  a  function 
of  signal  to  noise  ratio,  jamming  (or  interference)  power  to 
signal  power,  and  E,  which  is  defined  by  Eq.  (4.12).   As 
previously  stated,  E  is  the  ratio  of  interference  bandwidth 
to  receiver  bandwidth.   Figs.  4.15  and  4.16  show  the  performance 
of  the  white  noise  receiver  for  SNR  =  10  and  1,  respectively, 
at  various  values  of  JSR  and  E.   If  we  compare  these  figures 
to  Figs.  4.13,  4.14  which  show  the  performance  of  the  colored 
noise  receiver  that  has  been  optimized  to  the  specific  inter- 
ference, we  reach  the  conclusion  that  the  performance  of  both 
receivers  is  almost  the  same.   In  fact,  due  to  the  limited 
resolution  of  the  figures,  one  can  hardly  notice  any  differ- 
ence in  performance  at  all.   In  order  to  show  the  actual 
differences  in  performance,  Tables  4.1-4.5  present  numerical 
values  for  the  performance  of  both  receivers  under  various 
conditions.   The  tables  show  that  the  colored  noise  receiver 
always  has  better  performance.   However,  this  performance 
improvement  is  in  the  order  of  a  few  percent  in  the  best 
cases . 
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Table  4.1 

SNR  =-10.0 

JSR  = 

=  0.1 

E 

Pe  Colored 

P   White 
e 

Ratio  of  P 

0.1 

0.0043 

0.0043 

0.9997 

0.5 

0.0031 

0.0031 

0.9950 

1.0 

0.0021 

0.0021 

0.9859 

1.5 

0.0014 

0.0015 

0.9774 

2.0 

0.0010 

0.0010 

0.9712 

2.5 

0.0007 

0.0008 

0.9676 

3.0 

0.0006 

0.0006 

0.9661 

3.5 

0.0004 

0.0004 

0.9661 

4.0 

0.0003 

0.0003 

0.9671 

4.5 

0.0003 

0.0003 

0.9689 

5.0 

0.0002 

0.0002 

0.9709 

Table  4.2 

SNR  =10.0 

JSR  = 

=  1.0 

E 

Pe  Colored 

P   White 
e 

Ratio  of 

0.1 

0.1594 

0.1595 

0.9993 

0.5 

0.1434 

0.1450 

0.9892 

1.0 

0.1249 

0.1284 

0.9728 

1.5 

0.1089 

0.1136 

0.9586 

2.0 

0.0952 

0.1005 

0.9477 

2.5 

0.0835 

0.0889 

0.9400 

3.0 

0.0736 

0.0787 

0.9348 

3.5 

0.0651 

0.0698 

0.9315 

4.0 

0.0577 

0.0  621 

0.9298 

4.5 

0.0514 

0.0553 

0.9292 

5.0 

0.0459 

0.0493 

0.9294 
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SNR  =10.0 


0.1 
0.5 
1.0 
1.5 
2.0 
2.5 
3.0 
3.5 
4.0 
4.5 
5.0 


P   Colored 
e 


0.3733 
0.3631 
0.3511 
0.3401 
0.3298 
0.3203 
0.3114 
0.3030 
0.2951 
0.2876 
0.2805 


Table  4.3 

JSR  =10.0 

P   White 
e 

Ratio  of  P 

0.3736 

0.9991 

0.3657 

0.9930 

0.3560 

0.9865 

0.3465 

0.9816 

0.3372 

0.9780 

0.3283 

0.9755 

0.3198 

0.9738 

0.3115 

0.9726 

0.3036 

0.9719 

0.2960 

0.9715 

0.2887 

0.9713 

Table  4.4 

SNR  =1.0 

JSR  =  10.0 

E 

P   Colored 
e 

P   White 
e 

Ratio  of  P 

0.1 

0.3763 

0.3763 

0.9999 

0.5 

0.3681 

0.3690 

0.9978 

1.0 

0.3580 

0.3600 

0.9945 

1.5 

0.3484 

0.3513 

0.9917 

2.0 

0.3394 

0.3429 

0.9897 

2.5 

0.3311 

0.3350 

0.9884 

3.0 

0.3233 

0.3274 

0.9875 

3.5 

0.3161 

0.3202 

0.9871 

4.0 

0.3093 

0.3134 

0.9869 

4.5 

0.3029 

0.3069 

0.9870 

5.0 

0.2969 

0.3008 

0.9872 
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Table    4.5 


SNR   =1.0 


JSR   =1.0 


E 

P   Colored 
e 

P   White 
e 

Ratio  of  P 

0.1 

0.2033 

0.2033 

1.0000 

0.5 

0.1936 

0.1937 

0.9993 

1.0 

0.1826 

0.1830 

0.9979 

1.5 

0.1730 

0.1735 

0.9967 

2.0 

0.1647 

0.1653 

0.9959 

2.5 

0.1574 

0.1582 

0.9954 

3.0 

0.1512 

0.1519 

0.9952 

3.5 

0.1457 

0.1464 

0.9952 

4.0 

0.1410 

0.1416 

0.9954 

4.5 

0.1368 

0.1374 

0.9957 

5.0 

0.1330 

0.1336 

0.9960 
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One  must  conclude  that  under  the  given  conditions,  the  colored 
noise  receiver  does  not  perform  significantly  better  than  the 
white  noise  receiver.   The  next  logical  step  is  to  choose 
optimum  signals  and  determine  whether  the  colored  noise 
receiver  affords  greater  performance  improvement  over  a  white 
noise  receiver  operating  in  the  same  environment. 

G.   PERFORMANCE  OF  THE  "COLORED  NOISE  RECEIVER"  WITH 
OPTIMUM  WAVEFORMS 

As  has  been  demonstrated  the  colored  noise  receiver  does 
not  perform  significantly  better  than  a  white  noise  receiver 
for  the  signal  choice  of  the  previous  section.   However,  the 
performance  of  a  colored  noise  receiver,  unlike  that  of  a 
white  noise  receiver,  depends  on  the  signal  waveforms.   In 
this  section  we  use  the  optimum  signal  waveforms  derived  in 
Section  D  (Eq.  (4.20))  and  analyze  the  performance  of  the 
receiver  which  is  designed  to  match  these  waveforms. 

The  performance  of  the  colored  noise  receiver  with 
equiprobable  binary  antipodal  signals  is  given  by 


!    T  1/2 

Pe   =   ERFC*(i   /   yd(t)hd(t)dt)  (2.26) 


where  y^(t)  was  defined  by  Eq.  (4.20),  namely 


y,(t)   =   2A  sinbt  .  (4.20) 

a 
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The  correlating  signal  h,  (t)  is  defined  by  Eq.  (B.6)  and 
repeated  here  for  convenience 

hd(t)   =   C  sinbt  +  CK1£Yt  +  CK2e~yt  (B.6) 

where  the  constants  C,  K, ,  and  K2  are  defined  by  Eqs.  (4.22), 
(4.23),  (4.24)  respectively. 

Substituting  Eqs.  (4.20)  and  (B.6)  into  the  integral  of 
Eq.  (2.26) ,  yields 


P    =   E 

0 


1    T  +-         *-     X/2 

RFC*(±-   /  2A  sinbt  [C  sinbt  +CK1eyr  +  CK  e  Y  ]dt) 

(4.42) 


Evaluating   the   integral   yields 


Pe      =      ERFCjAC[j(l  -512bTbT)     +   -f^-jtY   sinbT-b  cosbT] 

Y   +  b 

K1be~YT     K2b     K2e~YT  1/2 

+   — 2 2  +  ~2 2  ~  ~~ 2 2^   sinbT  ~b  cosbT)  ]  } 

y   +b     y   +t)     Y   +b 


(4.43) 


We  define  a  normalized  frequency  b, 


b   =  £■  (4.44) 


The  normalized  frequency  b  is  the  ratio  of  the  signal  frequency 
to  the  noise  bandwidth.   Notice  that  when  dealing  with  bandpass 
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signals,  this  is  the  ratio  between  the  frequency  of  the 
modulation  to  the  noise  bandwidth. 

Using  Eqs .  (4.44)  and  (4.12),  the  factor  bT  can  be 
written  as 


bT   = 


T   =   bE  . 


(4.45) 


Also,  using  Eqs.  (A. 17)  and  (4.11)  yields 


yT   =   itu  6  T   =   m,E 


(4.46) 


Substituting   Eqs.     (4.46),     (4.45),     (4.44)    and    (4.22)     into 
Eq.     (4.43)    yields 


4A    T,,      sin  2bE,  ,,  ,rr2>    r- 
-d )  (1+b    ) 


P       =     EKPC., 


Nl 


2bE 


N0  -2 


-Em. 


1  +  2K, 


m,  sin  bE-b  cosbE  +  be 

~T  2  ,  r-2-,  M        sin  25E-, 
-     Elml+b  ][1  ~       2EE     ] 


1-, 


-Em, 


+  2K, 


b  -  e         (rru  sinbE  +b  cosbE) 
„ ,  2  ,  r-2>  ,,      sin  2bE, 

.  E(mi+b  )(1 — mr-] 


1/2 


i 


(4.47) 


The  signal  energy  E   is  given  by 


AT.     sin  2bE, 


2bE 


(4.48) 


Substituting  Eqs.  (4.48)  and  (4.40)  into  Eq.  (4.47)  yields 
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P       =     ERFC     '         2SNR 


e  *  ^  4  (JSR)  (SNR)      , 


(    E(l+b2) 


_  _  -Em.. 

m,  sinbE  -b  cosbE  +b£ 

1  +  2K,   — 

__, ,    2  ,  r-2.  ,,        sin  2bE» 
E(m,  +b  )  (1 = ) 

2bE 


_       _Emi  _  )    1/2 

b  -  £  (m,  sinbE  +b  cosbE)  ^ 

+  2K2- 


E(nL2  +b2)  (1  -  sinibE) 
2bE 


I 


(4.49) 


The  constants  K.  ,  K?  are  defined  by  Eqs.  (4.23)  and  (4.24). 
Eq.  (4.49)  is  a  rather  formidable  expression.   It  can  be 
calculated  by  a  computer. 

Figs.  4.17  and  4.18  show  the  performance  of  the  colored 
noise  receiver  matched  to  sinusoidal  waveforms,  for  SNR  =  10 
and  E  =  1.0  as  a  function  of  b.   They  also  show  for  comparison 
purposes  the  performance  of  the  colored  noise  receiver  matched 
to  rectangular  pulse  waveforms.   The  important  conclusion  one 
can  draw  is  that  by  increasing  b,  the  performance  of  the 
colored  noise  receiver  matched  to  sinusoidal  waveforms  improves 
significantly . 

Consider  the  following  numerical  example  by  referring  to 
Fig.  4.17.   For  SNR  =  10  and  JSR  =  1  the  colored  noise  receiver 
matched  to  rectangular  pulses  has  P   =0.15.   The  colored 

noise  receiver  matched  to  sinusoidal  pulses  has  a  P   =  0.065 

c  e 

for  the  same  JSR  and  SNR  values  and  b  =  1.   This  represents 

a  significant  improvement  however  neither  receiver  can  function 

properly  at  such  a  high  P  .   Increasing  the  signal  frequency 

_  -3 

up  to  b  =  6  causes  P   to  drop  to  a  value  of  10   .   Now  the 
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latter  receiver  can  operate  properly  even  under  severe 
jamming  conditions.   Notice  that  in  spite  of  the  fact  that 
the  signal  frequency  has  been  increased,  the  signal  energy 
collected  by  the  receiver  is  not  decreased  by  as  much  as  that 
of  the  interference.   Thus  receiver  performance  improves. 
One  must  conclude  that  it  is  best  to  increase  signal  frequency 
as  much  as  possible  until  hardware  constraints  are  reached. 
Figs.  4.19  and  4.20  show  data  similar  to  that  presented  in 
Fig.  4.17  except  that  now  the  interference  bandwidth  is  much 
smaller.   Nevertheless,  similar  conclusions  can  be  reached. 

H.   PERFORMANCE  OF  THE  WHITE  NOISE  RECEIVER  WITH  SINUSOIDAL 
PULSES 

As  was  demonstrated  in  the  previous  section,  the  use  of 
sinusoidal  optimum  pulses  significantly  improved  the  per- 
formance of  the  colored  noise  receiver  in  comparison  to  the 
performance  when  rectangular  suboptimum  pulses  are  used.   In 
this  section  we  evaluate  the  performance  of  the  white  noise 
receiver  with  sinusoidal  pulse  signals  in  order  to  determine 
whether  the  improvement  discussed  above  for  the  colored  noise 
receiver  also  occurs  for  the  white  noise  receiver.   The  white 
noise  receiver  performance  in  the  presence  of  white  noise  is 
independent  of  the  signal  waveform.   However  when  an  additional 
(colored)  interference  is  introduced,  the  performance  of  the 
receiver  is  affected  by  the  signal  waveforms  as  will  be 
demonstrated  in  the  sequel. 
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The  inputs  to  the  receiver  are: 

Hypothesis  HQ : "    z(t)   =   yQ(t)  +  nc(t)  +  n  (t) 

0  1  t  1  T 
Hypothesis  H^i  z(t)   =   y1(t)  +  n  (t)  +  n  (t) 

The  signal  waveforms  are 

y, (t)   =   A  sin  bt 


so  that 


yQ (t)   =   -A  sin  bt 


y,(t)   =   2A  sin  bt  (4.51 


(4.50) 


The  energy  per  bit  is  given  by 


T  T  T 

E        =        j      y?(t)dt   =    /  y2(t)dt      /   A2  sin2btdt 
s      0  0    u  0 


2 

A  Tr,    sin  2bT,  ,.    cnx 

"2-[1 2bT— ]  (4'52) 


The  statistic  £  generated  by  the  receiver  has  a  conditional 
mean  given  by  Eqs.  (4.29)  and  (4.30). 

The  conditional  variance  is  given  by  substituting  Eqs. 
(4.33),  (4.34),  (4.35)  and  (4.50)  into  (4.32).   Performing  the 
substitution  yields 
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mm  MR  NT 

Var[£|H    ]       =      =jL     /        /      A2  sinbt  sinbx  [-2-  e"6  ' t_T  '    +  -i  6(  t-T  )  ]  dtdr 
1  s    0        0  q  Z 


Nl  T         T  N    B  rI*    T| 

-±  +      /         /A     sinbt  sinbx    -^-  e    P  '  '  dtdx 

^00  4 


N  A2N    6  T  t  , 

-^  +      A1?u        /      sin  bt[    /      sin  bx    e"eit   T  '  dx 
^  s      0  0 


T 
+      /      sinbx    e"3 (x_t)dx]dt    .  (4.53) 

t 


Performing  the  integration  and  substituting  Eq.  (4.52)  yields 


VarU|H  ]   =   4  +  4  23   2 [1  +  ^silfbT  ]  (4'54) 

3   +bZ      E(l  -   %b%  ) 


where  E  is  defined  by  Eq.  (4.12). 

Applying  the  same  procedure  for  Var[£|Hn]  yields  the 
same  result  as  Eq.  (4.54) .   Thus 

Var[£|H1]   =   Var[£|HQ]  .  (4.38) 

The  performance  of  the  receiver  can  now  be  calculated 
from 
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=   ERFC. 


**S /         (4.56) 


e  *  <  /  N,   N_   D2  .  2,  _ 

1  +  0  _3.u  +     Bin  bT 


=   ERFC* 


2    2  $2+b2'    E(1.siybT, 


2Es/Nl 

1  +^0    bL(1  +    sin2bT 


N,  D2^,2  l     ^  ,,   sin  2bT, 
1  3  +b       E(l 2b^' 


Substituting  Eq.  (4.49)  yields  (4.57) 


P    =   ERFC.  2(SNR)  <  (4'57) 


i  ,  4(JSR)  (SNR)     1    n      sin2bT    , 

E  b  2 

1  +  (s-)  _M   sin_2bT 

3  E(1  "  2bT  ) 


Eq.  (4.57)  specifies  the  performance  of  the  white  noise  re- 
ceiver with  sinusoidal  pulses  in  the  presence  of  white  and 
colored  noise  interference,  as  a  function  of  SNR,  JSR,  E 

(defined  by  Eq.  (4.12)),  and  (q-)  which  is  the  ratio  of  the 

p 

signal  modulation  frequency  to  the  bandwidth  of  the  inter- 
ference.  Figures.  4.21-4.25  show  a  comparison  of  the  perfor- 
mance of  the  receivers  analyzed.   That  is,  the  white  noise 
receiver  and  the  colored  noise  receiver,  both  with  sinusoidal 
pulse  waveform  and  with  rectangular  pulse  waveform  inputs. 

By  analyzing  these  figures,  one  can  reach  the  following 
conclusions : 

1.   For  narrow-band  interference  like  the  one  shown  in 
Fig.  4.21  (E  =  0.1),  the  colored  noise  receiver  matched  to 
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sinusoidal  pulses  performs  much  better  than  the  white  noise 
receiver  with  a  similar  input.  Let's  define  an  improvement 
factor,  I, 


I   = 


Pe,2 


where  PQ  -,  is  the  performance  of  the  former  receiver  and  P 

e,l        c  e, 2 

is  the  performance  of  the  latter  receiver. 

Looking  at  Fig.  4.21,  one  can  see  that  for  low  modulation 
frequency  (b) ,  the  improvement  factor  is  significant  and  can 
reach  a  value  of  30.   As  the  modulation  frequency  increases, 
the  improvement  factor  decreases.   At  very  high  modulation 
frequencies  both  receivers  have  almost  equal  performance  as 
can  be  seen  in  Figures  4.23,  4.24,  4.25. 

2.  For  low  modulation  frequencies,  the  white  noise  receiver 
with  sinusoidal  pulse  input  performs  worse  than  the  white 
noise  receiver  with  rectangular  pulse  inputs. 

3.  When  JSR  increases  without  bound,  both  receivers  are 
driven  into  saturation  and  the  improvement  factor  decreases 
as  can  be  seen  from  Fig.  4 .  22  . 

4.  As  the  bandwidth  of  the  interference  increases,  the 
improvement  factor  decreases  as  can  be  seen  from  Figures  4.21, 
4.23,  4.25.   However  the  colored  noise  receiver  exhibits 
better  performance  than  the  white  noise  receiver  with  similar 
inputs . 
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V.   RF  PREFILTER — COLORED  NOISE  THEORY  ANALYSIS 

A.  INTRODUCTION 

As  discussed  in  Chapter  II,  a  common  source  of  nonwhite 
Gaussian  noise  in  the  communication  channel  is  the  presence 
of  a  bandpass  element  between  the  transmitter  and  the  signal 
processing  sections  of  the  receiver.   The  most  common  such 
bandpass  element  is  a  low  noise  RF  preamplifier  used  to 
improve  the  sensitivity  of  the  receiver. 

In  this  chapter  we  analyze  the  effect  of  this  RF  preampli- 
fier, using  some  of  the  methods  and  results  previously  derived 
The  analysis  will  be  done  for  two  different  cases. 

1.  The  ideal  case  in  which  the  noise  figure  of  all  the 
elements  in  the  receiver  is  equal  to  unity. 

2.  The  more  realistic  case  in  which  the  receiver  components 
have  noise  figures  that  are  greater  than  unity. 

B.  THE  MODEL 

In  this  chapter,  two  receivers  will  be  analyzed  and  their 
performances  will  be  compared.   The  first  receiver  is  a 
binary  coherent  digital  receiver.   This  receiver  is  optimum 
for  discriminating  between  signals  received  in  an  additive 
white  Gaussian  noise  environment.   The  input  to  this  receiver 
consists  of  the  information  signals  with  the  additive  white 
Gaussian  noise  having  P.S.D.  level  N-/2.   This  receiver  is 
described  by  Fig.  5.1.   The  second  receiver  is  described  by 
Fig.  5.2.   It  consists  of  an  RF  preamplifier  at  the  front  end 
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White  noise  receiver 


Signal 

yx(t>- 

or 

y0(t) 


n  (t) 
w 


Figure  5.1 


Colored  noise  receiver 


Figure  5.2   Preamplifier  Receiver 
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of  the  receiver.   We  assume  that  the  preamplifier  has  a 
transfer  function  given  by 


H(f)  ■  W^fre  f5-1' 


The  input  noise  power  spectral  density  level  is  Nn/2.   The 
RF  preamplifier  produces  a  colored  noise  component  due  to 
the  white  input  noise,  having  power  spectral  density  given  by 


N        2 

4>"c(s)   =  -j   G2  — | j    (s  =  J27Tf)         (5.2) 

-s  +  3 


It  is  reasonable  to  assume  that  at  the  input  of  the  correlator 
there  is  also  an  additive  white  Gaussian  noise  component  due 
to  the  front  end  thermal  noise.'  The  total  noise  P.S.D.  at 
the  input  to  the  correlator  is  thus 


N     N         2 

<j>(s)   =  $  (s)  +  »  (s)   =  ~y   +  ~y   G2  \ j         (5-3) 

-s  +3 


Observe  that  the  information  signals  are  distorted  by  the 
preamplifier  and  therefore  the  output  of  the  preamplifier  no 
longer  produces  signals  y.(t),or  y, (t) ,  but  rather  yA(t)  or  y|(t) 
Clearly  y-l(t)  and  yA(t)  are  the  result  of  convolving  y,  (t) 
and  yn(t)  respectively  with  the  preamplifier  impulse  response. 

The  receiver  described  in  Fig.  5.2  is  optimum  for  dis- 
criminating between  yA(t)  and  y' (t)  provided  h, (t)  is  appro- 
priately chosen.   The  autocorrelation  function  of  the  total 
noise  is  given  by 
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N0         N0    2 
K(t)   =  -f  6  (t)  +  -%■   3GZ  exp(-e|x|) 


The  power  of  the  colored  noise  component  is  given  by 

N0    2 


Observe  that  unlike  the  model  of  Chapter  IV,  P   is  no  longer 
constant.   When  3  changes,  P   is  changed  also. 

C.   RECTANGULAR  PULSE  RECEIVER  DESIGN 

Assume  the  input  signals  to  the  RF  preamplifier  to  be 


0  <_   t  <  T 
Y1(t)  (5.5) 

T  <  t,  t  <  0 


to 


/  -A     0  <  t  <  T 
t)   = 


y0(t)   =  (5.6) 


! 


0      T  <  t,  t  <  0 


so  that 


/  2A     0  <_   t  <_   T 

yd(t)  =  (5.7) 

(  0      T  <  t,  t  <  0 


Since  the  impulse  response  of  the  preamplifier  is 


h(t)   =   Ge  3tu(t) 


91 


then 


2AG(1  ~e"et)     0  <_   t  <    T 

y^(t)  (5.8) 

(  0  t  <  0,  t  >  T 


where 


y^(t)   =   h(t)  *  yd(t) 

The  correlating  signal  in  the  colored  noise  receiver  is  the 
solution  to  the  integral  equation 


-j   hjt)  +   /   B  -£  G  exp(-3|t-u|)hd(u)du  =  y'  (t)  =  2AG(l-e3t) 
0 

(5.9) 


The  solution  to  this  Fredholm  equation  is  somewhat  more  compli- 
cated than  the  one  worked  out  in  Appendix  A  due  to  the  fact 
that  the  function  on  the  right  hand  side  of  the  integral 
equation  is  no  longer  a  constant  whose  derivatives  are  zero, 
but  rather  an  exponential  function.   The  detailed  solution 
to  Eq.  (5.9)  is  worked  out  in  Appendix  E. 

The  correlating  signal  h,  (t)  is  given  by 

+  Bm,  t       -3m,  t 

h,(t)   =  C  +  CK,£      +  CK„e  (5.10) 

a  1  z 


where 
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C   =      4A%  (5.11) 

NQ(1  +G  ) 


v^T? 


ml   =        G  (5.12) 


-Em,    in,  -1   -2Em, 

(ml+1) [£     "  iTTT  £     ] 

Kx   = (5.13) 

m,+l   m,-l   -2Em, 

m,  -1   m, +1 


-Em.,    m.+l 

(m-,-1)  [e    l    -  -A-T-] 

1  m1  -1 

K2   =     m~fl    ET^l   -2Em1  (5.14) 

m,  -1   m, +1 


Ignoring  the  constant  of  proportionality,  it  is  easy  to  see 
that  h,(t)  is  a  function  of  G  and  E  only.   Recall  that  G  is 
the  gain  of  the  preamplifier  and  E  is  the  ratio  of  the  preampli- 
fier bandwidth  to  the  receiver  bandwidth. 

Fig.  5.3  shows  a  plot  of  h,  (t)  versus  time  normalized 
to  the  pulse  width,  for  various  values  of  E  and  G  =  20  db. 
In  practical  design,  however,  the  preamplifier  bandwidth  will 
not  be  much  different  than  the  correlator  bandwidth.   Fig.  5.4 
shows  again  h,  (t)  for  various  values  of  E  with  G  =  0  db.   The 
white  noise  component  in  this  case  is  dominant  causing  the 
correlating  signal  h ,  (t)  to  be  almost  flat.   This  is  in  agree- 
ment with  known  results  on  the  correlation  operation  for 
receivers  operating  in  white  noise  interference  only  with 
constant  pulse  input. 
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D.   RECEIVER  PERFORMANCES — THE  IDEAL  CASE 

In  this  section  the  performance  of  the  receivers  shown 
in  Figs.  5.1  and  5.2  is  analyzed. 

The  performance  of  the  colored  noise  receiver  with  equi- 
probable  binary  signals  is  given  by  Eq.  (2.26). 


l    T 
=   ERFC*(±-   /   y'  (t)h,  (t)dt) 


1/2 


dv  -   d 


(2.26) 


If  we  substitute  y^(t)  as  given  by  Eq.  (5.8)  and  hd(t)  as 
given  by  Eqs.  (5 . 10 ) -  (5 . 14 )  ,  we  get 


Pe   =   ERFC*<*  /     8A  G  2 

0    NQ(1  +GZ) 


a.  3m, t     -3m..  t    1/2 

(1  -  e'pr)  (1  +K1e  x  +K2e   l    )dt) 


(5.15) 


Evaluating  the  integral  yields 


P       =     ERFC 


9a2   2  .      -E  ™1E  ,  (VUE.  _      "mlE 

2A  G  T     ..     1-e  e     '-l     „     e  -1  .  „    1-e 

= — \L = his., 

NQ(1+GZ)  * 


1       n^E    "Kl     Edn^l)       +K2       En^ 

ll/2 


+  K 


-E(m,+1) 
2       EOyKL)      ; 


(5.16) 


Since  A   T/N_    can   be    interpreted    as    the   SNR,    we   obtain 


P.     =     ERFC4 


2  -E  mlE 

2SNRG  1-c  b  e       -1 

E      +*1     it^E  1     E(m1-1)  *2       Ett^ 


(m,-l)E  -m,E 

e  -1  .  „     1-e 


(1-KT) 


-E(m,+1) 
*2       E(m  +1)      ; 


l/2 


5.17) 
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where  K,  and  K  are  defined  by  Eqs.  (5.13)  and  (5.14).  The 
white  noise  receiver  performance  for  equiprobable  antipodal 
signals  is  given  by 


P    =   ERFC* [/2SNR] 


Fig.  5.5  shows  a  performance  comparison  between  the  two  re- 
ceivers for  G  =  50.   It  is  clear  that  both  receivers  perform 
equally  as  well.   The  addition  of  the  preamplifier  did  not 
improve  the  performance  of  the  receiver  of  Fig.  5.2. 

E.   RECEIVER  PERFORMANCES — PRACTICAL  CASE 

In  the  previous  section  we  assumed  the  system  analyzed 
consisted  of  only  ideal  components.   These  components  didn't 
contribute  any  noise  to  the  overall  system. 

In  this  section  we  discard  this  assumption  and  instead 
work  with  practical  elements,  so  that  the  white  noise  receiver 
has  an  input  noise  figure  NF, .   We  may  thus  state  that  the 

power  spectral  density  level  due  to  the  input  white  Gaussian 

N0 
noise  is  no  longer  N-/2  but  rather  -=-  •  NF,  .   Therefore,  the 

performance  of  the  white  noise  receiver  for  equiprobable 

antipodal  signals  is 


P    =   ERFC* [/2SNR/NFT] 


The  preamp  receiver  has  also  a  non-unity  noise  figure  asso- 
ciated with  it,  which  we  denote  NF? .   Since  the  receiver  input 
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stage  is  assumed  to  be  a  low  noise  preamplifier  (otherwise 
there  is  no  benefit  in  utilizing  this  preamplifier) ,  we  can 
assume  that 


NF„   <   NF,  . 


The  performance  of  the  colored  noise  receiver  is  no  longer 
given  by  Eq.  (5.17)  but  must  be  modified  accordingly  to  yield 


P   =  ERFCg 


2SNRG 


NQ(NF,  +NF2G  ) 


,    -e        "iE  ,      (V1): 


-1 


m,E 


(n^-1) 


-m,E 


•m. 


-E 


+  K2(i% 


Oy-1) 


+ 


E(m,+1) 


)] 


1/2 


(5.18) 


Fig.  5.6  shows  the  performance  of  the  two  receivers  analyzed 
in  this  section.   It  is  clear  that  utilization  of  a  low  noise 
preamplifier  improved  significantly  the  performance  of  the 
colored  noise  receiver.   However  it  must  be  pointed  out  that 
if  we  design  a  receiver  which  consists  of  white  noise  re- 
ceiver and  a  preamplifier,  .ts  performance  would  almost  be  the 
same  as  the  performance  of  the  colored  noise  correlator  with 
preamplifier  as  described  in  Fig.  5.6.   The  major  contribution 
to  improved  performance  is  due  to  the  low  noise  amplifier  and 
not  due  to  the  specific  correlator  used  in  the  system. 
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VI.   JAMMING  THE  COLORED  NOISE  RECEIVER 

A.  INTRODUCTION 

In  previous  chapters  the  design  and  performance  of  re- 
ceivers operating  in  an  environment  consisting  of  colored 
noise  interference  were  analyzed.   The  interference  was 
assumed  known  and  the  receiver  was  optimized  to  the  presence 
of  that  interference.   In  this  chapter  the  model  is  expanded 
by  considering  the  presence  of  a  hostile  jammer  attempting  to 
jam  the  communication  channel.   The  jammer  is  hostile  in  the 
sense  that  its  parameters  are  not  known  to  the  receiver  designer 
The  main  purpose  here  is  to  find  the  optimal  jammer  waveform 
that  causes  maximum  damage  to  a  communication  channel  of  the 
type  analyzed  in  previous  chapters  and  to  determine  whether 
the  colored  noise  receiver  is  more  or  less  sensitive  than 
a  white  noise  receiver  to  jamming  signals. 

B.  THE  MODEL 

The  system  model  is  described  in  Fig.  6.1.   It  consists 
of  a  digital  coherent  communication  receiver  operating  in 
the  presence  of: 

1.  Colored  noise  interference 

2.  Additive  white  Gaussian  noise 

3.  Jamming  signal. 

The  digital  information  is  transmitted  via  binary,  baseband 
signals  that  encounter  baseband  interference  and  jamming. 
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Extensions  to  bandpass  signal  analysis  is  straightforward 
in  light  of  the  results  described  in  Chapter  III,  Section  C. 

The  P.S.D.  of  the  colored  noise  interference  is  assumed 
to  be  the  same  as  in  Chapter  IV.   Specifically,  the  colored  noise 
P.S.D.  used  here  is  given  by  Eq.  (4.2).   Also  the  signals  y, (t) 
and  yQ(t)  are  defined  by  Eqs .  (4.8)  and  (4.9).   The  jammer  is 
assumed  to  be  deterministic  and  its  specific  waveshape  will 
be  determined  in  the  solution  to  the  optimization  problem. 

C.   DERIVATION  OF  THE  OPTIMAL  JAMMING  WAVEFORM 

The  decision  process  due  to  the  presence  of  a  jammer 
becomes 


Hypothesis  Hn :   z(t)   =   yA (t)  +  n  (t)  +  n  (t)  +  n.(t)     (6.1) 

0  0        C        w        J 

0  £  t  <_   T 

Hypothesis  Hn :   z(t)   =   y,  (t)  +  n  (t)  +  n  (t)  +  n.(t)     (6.2) 

1  1       c       w       j 


Observe  that  n.(t)  is  modeled  here  as  a  deterministic  waveform. 
The  signals  yn (t)  and  y, (t)  are  assumed  to  be  antipodal.   The 
energies  of  the  signals  are  given  by 


T  T 

Ec  =   /  y?(t)dt  =   /  y^(t)dt  .  (6.3) 

s     0  0    u 


The  receiver  generates  the  statistic  I   which  is  a  Gaussian 
random  variable.   Its  conditional  mean  is  given  by 
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U|H    ]       =      E{    /      h     (t)  [y(t)    +n    (t)    +  n.(t)    +n(t)]dt} 

J-  i-v  Q.  J_  C  J  W 


/      h    (t)  [y    (t)    +   n.  (t)]dt      =      mn 


(6.4) 


and   similarly 


EU|HQ]       =         /      hd(t)[yQ(t)    +   n.(t)]dt      = 


A 

=      m. 


(6.5) 


The  conditional  variance  of  I    is  given  by 


VarUl^} 


Var{£|HQ} 


=   E 


h,  (t)  [n  (t)+n  (t) ]dt 
d     w     c 


A   ^2 


(6.6) 


Note  that  since  the  jammer  waveform  is  modeled  as  determinis- 
tic, it  does  not  affect  the  variance  of  the  statistic  I. 
This  would  not  be  the  case  if  the  jammer  waveform  were  a 
sample  function  of  a  random  process.   The  receiver  threshold 
Y  for  antipodal  equiprobable  signals  is  equal  to  zero  as 
previously  demonstrated  (see  Eq.  (2.1)). 
Define  now 


£  -  y  /  y.(t)h  (t)dt 


(6.7) 
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2 
From  Eq.  (2.23),  we  know  that  2E,   =   oQ.      Introducing  the  short- 
hand notation 


T 
(h  ,n  )   =    /   h  (t)n  (t)dt  (6.8) 

a   -1       0    a     J 


the  performance  of  the  receiver  of  Fig.  6.1  can  be  shown 
to  be  given  by 


2  [-?-(hd'ni)]/a£  2 

e  z   [?-(hd/n.)]/a£      /2?  *   -oo  /2?  dx 


(6.9) 


,               £-(h,,n.)         ,               5+(h,,n.) 
=     i-  ERFC .  [ — J_]    +  i.  erfc    [ E_J_] 

2  /2T  2  /2T 


where  P(Hn)   =   P{H, }   =   y  nas  heen  assumed. 
Define 


a   =   (h,,n.) 

J 


and  differentiate  P   with  respect  to  a.   This  yields 


-^    =    £ p =-£ ]     =    - sinh/C72a         (6.11) 

9a  /27F2T  /4^ 
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Since  £  is  non-negative,  it  is  clear  that 


>  0    a  >  0 


_E  =     =o    a  =  0  (6.11) 

3a 

<  0    a  <  0 


Therefore  a  =  0  is  a  minimum  point  of  P  ,  and  by  making 

[  ot  1  as  large  as  possible,  P   is  maximized,  because  Eq.  (6.11) 

shows  that  P   is  monotonic  in  |a|. 

From  the  Cauchy-Schwarz  inequality,  it  can  be  seen 
that 


a      =      (hd,n.)   <   ||hd||  •  Iln_.ll  (6.12) 

with  equality  holding  if 


n.  (t)   =   K  h,(t)  (6.13) 


where  K  is  an  arbitrary  constant. 

The  energy  E    of  the  jammer  is  given  by  (from  Eq.  (6.13)) 

n] 


T 

E  .   =   K2  (      h2(t)dt   =   K2.||h,||2  (6.14) 

nj        J         d  ' '  d1 ' 


Therefore,  we  must  have 


/E  . 
K   =   — 21 (6.15) 

MhdH 
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and  the  optimum  jammer  that  is  energy  constrained  is  given 
by 


/IT 
n  (t)   =   ^_  h,  (t)  (6.16) 

:       llhdH 

This  derivation  is  valid  for  both  a  white  noise  and  a  colored 
noise  coherent  receiver.   The  only  difference  is  in  h, (t) . 
For  a  white  noise  receiver,  h , (t)  is  proportional  to  yH(t). 
Therefore  the  waveform  of  the  optimum  jammer  will  be  related 
to  the  waveforms  used  to  transmit  the  binary  information. 
However  for  the  colored  noise  receiver,  h, (t)  is  no  longer 
directly  related  to  the  signals  y, (t)  and  yn(t).   Thus  the 
optimum  jammer  may  have  a  waveshape  that  has  no  resemblance 
to  the  waveforms  used  to  transmit  the  binary  information.   It 
is  clearly  feasible  to  implement  an  optimum  jammer  against 
a  white  noise  receiver.   All  that  must  be  done  is  to  transmit 
the  difference  of  the  signal  waveforms,  or  use  a  repeater 
channel  [Ref.  15] .   However  it  is  almost  impossible  to  opti- 
mize a  jammer  against  a  colored  noise  receiver  unless  all  the 
details  about  the  correlator  receiver  being  used  are  known. 

D.   PERFORMANCE  OF  THE  WHITE  NOISE  RECEIVER  WITH  OPTIMAL 
JAMMING 

In  order  to  properly  evaluate  the  effect  of  jamming  on 

colored  noise  receivers,  it  is  necessary  to  first  evaluate 

the  effect  of  jamming  on  white  noise  receivers.   The  results 

on  the  latter  can  then  be  used  as  a  reference,  to  which  results 
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on  the  former  can  be  compared.   We  assume  here  the  model 
described  in  Fig.  6.1  with  the  only  difference  being  that 
now  a  white  noise  receiver  is  used  in  place  of  the  colored 
noise  receiver.   That  is,  h,(t)  must  be  replaced  by  y,(t). 
The  receiver  performance  without  jamming  (however  with  noise 
interference)  has  been  evaluated  in  Chapter  IV,  Section  F. 
A  statistic  I    is  generated  by  the  receiver  of  Fig.  6.1,  where 
the  conditional  moments  of  the  statistic  are  given  by  Eqs. 
(4.29),  (4.30),  (4.31)  and  repeated  here  for  convenience. 

EU|H,  ]   =   /E~  (4.29) 

-L  o 


EU|H0]       =      -    /Es  (4.30) 

VarU|Hl]       =      VarU|H0]       =      ^  +?^(E-l+exp(-E)  } 

"      °l  (4.31) 


Observe  that  E  is  defined  by  Eq.  (4.12). 

The  performance  of  this  receiver  is  given  by  Eq.  (4.39), 
namely 


E  /E 

-^ }   =   ERFC*  {/-§-}       (4.39) 

'l  |  IN0[E-l+exp(-E)]  V  a£ 


for  equiprobable  signals.   Here  again  the  deterministic  jammer 
affects  the  conditional  means  of  1   but  not  its  conditional 
variances. 
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We  now  obtain 


T 
EU|H  -]   =    /   h  (t)  [y  (t)  +n.(t)  ]  dt 
0  J 


T  T 

/   h  (t)y,(t)dt+   /   h,(t)n.(t)dt  (6.17) 

0    d    L  o    d    : 


T  T 

E[£|H  ]   =   -   /   h  (t)y  (t)dt  +   /   h  (t)n.(t)dt  (6.18) 

0  0  D 


Substituting  the  optimum  jammer,  derived  in  Section  C  and 
given  by  Eq.  (6.16),  yields 


T  y,(t)               T        /E  .  y,(t) 

E[i\E       ]     =      /     _S_y(t)dt+     /     y   (t)    .."J.  |   -5 dt 

1       0    2/E"    X  o     d     I  ^d1  I  2/E~ 

s  s 


=   /E   +  /E  .  (6.19) 

s     nj 


E[£|Hn]   =   /E   -  /E  . 
1  0        s     nj 


Recall  that  for  the  white  noise  receiver  h ,  (t)  =  y,(t)/2/E 

d      2  d       s 

The  performance  of  the  white  noise  receiver  becomes 


,                     /E                  /E     .               ,                     /E  /E     . 

P    =      ±-  ERFC*[ -(1    +    — ^l)]  +   =■  ERFCA[— ^(l ^-)  ]  (6.20) 

8    2         a£       /W                 2                      °l  /E~ 

s  s 


Observe  that  if  no  jammer  is  present,  Eq.  (6.20)  takes  on 

the  form  of  the  P   for  a  white  noise  receiver  in  the  presence 

e 
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of  noise  interference  only.   The  jammer  effect  can  be  seen  from 

the  introduction  of  the  (1  ± -)  factors. 

s 

E.   PERFORMANCE  OF  THE  COLORED  NOISE  RECEIVER  IN  THE 
PRESENCE  OF  JAMMING 

The  performance  of  the  colored  noise  receiver  in  the 

presence  of  jamming  was  derived  in  Section  C  and  the  receiver 

P   is  given  by  Eq.  (6.9),  namely 
e  "  -* 


!       5+<hd,n.)     1  S-(hd,n  ) 

P^   =   y  ERFCA[ - -L-]  +  ±-ERFC*[ ^-]  (6.9) 

e      2  /2T        *  /2f 


We  assume  that  the  model  described  in  Section  B  is  valid  in 
the  foregoing  analysis. 

The  performance  of  the  receiver  in  Fig.  6.1  when  no  jammer  is 
present  was  analyzed  in  Chapter  IV  and  its  P   given  by  Eq. 
(4.26) ,  namely 

-Em,   -2E(m,+l) 

2(m+l)-4e   X-2  e        (m -1)  1/; 

P   =  ERFC*(/£72)   =  EBFCA =^ T(l  + -, - =——  -H 

e        *   '  *  i,,T^/^^4         m,+l   itL-1  -2Em 


1+(JSR)  (SNR)i       _   .V1   V1  _^\ 


E 


2SNR  V2 


=  ERFC*[ ^^ ^f-^E,!^)]  (6.21) 


1+(JSR)  (SNR)- 


where  E  and  m,  were  defined  in  Eq.  (4.12)  and  Eq.  (4.17)  respectively 
and  JSR  is  the  ratio  of  interference  power  to  signal  power.   In 
order  to  evaluate  the  performance  of  the  receiver  analyzed 
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above  under  jamming  conditions,  the  factor  (h,  ,n.)  must  be 
evaluated  for  various  jamming  waveforms.   The  analysis  will 
be  carried  out  for  two  different  cases. 

1.  The  jammer  has  the  same  waveshape  as  the  information 
signal  waveforms  (i.e.,  rectangular  pulses). 

2.  The  jammer  has  been  optimized  according  to  the  results 
of  Section  D. 

Case  1: 

The  jammer  waveform  is  given  by 


n  (t)   =   —^1  y1it)  (6.22) 


/E 
s 


Substituting  Eq.  (6.22)  into  Eq.  (6.8)  yields 


T  T  /E  . 

(h,,n.)   =   /   h,(t)n.  (t)dt  =    /   — ^  y,  (t)  h ,  (t)dt 

d   3       o     d     :         0    /E"   X     d 

s 


/eT   T 
£1     j      y1(t)hd(t)dt       (6.23) 


'E    0 
s 


Substituting  Eq.  (6.7)  into  Eq.  (6.23)  yields 


/E  • 
(h,,n.)   =   —21  i  (6.24) 

d  3  /E~ 

s 


Furthermore,  substitution  of  Eq.  (6.21)  and  Eq.  (6.24) 
into  Eq.  (6.9)  yields  finally 
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Pe      =      i  ERFC^[/C/2(1    +    /Enj/Es)]     +   |   ERFC^[/?72(1    -    ^Enj/Eg) ] 


1    „M„     r      /  2SNR 


rERFC,[     /  —    -4    f^^)  (1    +    ^-TfT)] 


1   +  (JSR)  (SNR)^ 


x  pSNRfj^CE^    ) 


+    2    ERFC*U/ ± — Z-d    "    /En-i/E S)]  (6'25) 

'  1  +  (JSR)  (SNR)|        nj       S 


Analysis  and  simulation  carried  out  on  Eq.  (6.25)  reveals 
that  this  receiver  performs  at  almost  the  same  level  as  the 
white  noise  receiver  analyzed  in  Chapter  IV,  Section  F.   The 
performance  of  the  white  noise  receiver  without  jammer  is 
given  by  Eq.  (4.39).   Its  performance  is  almost  the  same  as 
that  of  the  colored  noise  receiver  without  jammer  whose  per- 
formance is  given  by  Eq.  (6.21).   This  fact  was  demonstrated 
by  the  numerical  results  presented  in  Tables  4.1-4.5. 

When  the  jammer  is  introduced,  the  arguments  of  the  error 
function  for  both  receivers  has  to  be  modified  by  the  same 


factor  (1  ± /E  .//e    ) .      Therefore  the  performance  of  the  two 
ny       s  ^ 

receivers  under  jamming  conditions  remains  almost  identical. 
Case  2; 

The  optimum  waveform  jammer  analyzed  in  Section  D  is 
given  by  Eq.  (6.16).   Substituting  Eq.  (6.16)  into  Eq.  (6.8) 
yields 


T   /eT  

(h,,n.)   =    /   2i_  h,  (t)h,  (t)dt  =  /E~  |  |h,  |  |  .      (6.26) 

d   D       n    I  Iv,  I  I   d    d  n:     d 


0   M  hd 
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Substituting  now  the  appropriate  h,  (t)  which  is  given  by 
Eqs.  (4.13)-(4.16)  yields 


-2  Em. 


^E        -ioxtt-,  <->  -Em,     -,  ,        1 

<W  =  ^  -jjJniMi'1  +  2K2£      +  K2(i^ — > 

v  t,         i  i — 1 

S  E 

(i-rErai)  ^2 

+  4K2  j — k—L]  (6-27> 


where  K2  is  defined  by  Eq.  (4.15). 

Equation  (6.27)  can  be  written  in  the  form 


/E- 


(Vnj>   "  ^   .     ,4(JSR)(SW  ^2^1^)  (6.28) 


S 


Substituting  Eq.s  (6.28),  (6.21)  into  Eq.  (6.9),  yields  the 
performance  of  the  colored  noise  receiver  in  the  presence  of 
the  optimum  jammer.   This  result  is 


,  /   2SNRf,  (Ei  ,m 7)  /E~7  /fTT13n7) 

t,  Sill 


,         I  2SNRf.(E1,mJ      /f2^1'ml) 

+  2  ERFC<  -  ,4(JSRnsNR)-(1  +  /V^Vvv7 )]    <6-29) 


where  f,  and  f~  are  defined  by  Eqs.  (6.21)  and  (6.28) 
respectively. 

In  Figures  6.2-6.5,  a  comparison  between  the  performance 
of  the  colored  noise  receivers  analyzed  under  the  two  jamming 
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conditions  is  shown.   The  pertinent  equations  are  Eq.  (6.29) 
and  Eq.  (6.25).   All  results  are  a  function  of  the  jammer 
to  signal  energy  ratio,  E  ./E    ,  denoted  as  JSR*.   The  figures 
show  that  the  optimum  jammer  (n  .(t)    =  Kh,(t))  causes  much  more 
damage  to  the  receiver  performance  than  the  suboptimum  jammer 
(n  .(t)   =   Ky,  (t)).   If  one  tries  to  determine  what  JSR  value  is 

J 

required  to  cause  the  receiver  to  operate  at  a  certain  P  , 
it  can  be  seen  that  when  the  optimal  jammer  is  used,  less 
jamming  power  is  needed  (somewhere  between  6-12  db  less 
jammer  power)  than  if  the  suboptimum  jammer  of  Case  1  is 
used. 
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VII.   CONCLUSION 

This  thesis  presents  several  applications  of  the  theory 
of  signal  detection  in  the  presence  of  colored  noise.   In 
it,  the  analysis  of  practical  design  implications  of  the 
theory  is  carried  out  and  the  evaluation  of  the  performance 
of  receivers  designed  according  to  this  theory  is  undertaken. 
The  design  of  a  digital  receiver  operating  in  a  colored  noise 
environment  requires  the  solution  of  a  specific  Fredholm 
integral  equation.   In  order  to  solve  the  Fredholm  integral 
equation,  the  designer  must  have  available  an  analytical 
expression  for  the  autocorrelation  function  of  the  colored 
noise  and  also  know  the  signal  waveforms  being  used  to  trans- 
mit the  digital  information.   Once  this  information  is  avail- 
able, solution  of  the  Fredholm  integral  equation  must  be 
attempted.   As  discussed  in  Chapter  III,  analytical  solutions 
do  not  always  exist.   Even  if  solutions  do  exist,  the  proce- 
dure for  solving  the  integral  equation  is  at  best  tedious. 
In  most  of  the  work  undertaken,  a  relatively  simple  case  in 
which  the  colored  noise  was  the  output  of  an  amplifier  stage 
followed  by  a  first  order  Butterworth  filter  was  analyzed. 
This  leads  to  an  analytical  expression  for  the  noise  autocorre- 
lation for  which  the  solution  of  the  integral  equation  exists 
and  is  tractable.   Several  cases  were  presented  in  Appendices 
A-E. 
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The  receiver  structure  and  correlator  signal  waveform 
were  derived  in  Chapter  IV  for  the  cases  in  which  the  binary 
information  was  transmitted  using  either  rectangular  or  sinu- 
soidal pulses.   As  shown  in  Appendix  D,  the  design  of  such  a 
receiver  was  quite  feasible  once  the  noise  autocorrelation 
was  specified  and  an  analytical  solution  to  the  Fredholm 
integral  equation  was  obtained.   Next/  the  performance  of  this 
receiver  was  evaluated  and  compared  to  that  of  a  white  noise 
receiver  operating  in  the  same  environment.   The  results 
showed  that  both  receivers  performed  almost  identically  with 
only  a  few  percent  difference  in  receiver  error  probabilities 
when  rectangular  pulses  were  used  for  signaling.   Only  when 
optimum  sinusoidal  pulses  were  used  in  place  of  the  rectangular 
pulses,  a  major  improvement  in  the  performance  of  the  colored 
noise  receiver  in  comparison  to  the  performance  of  the  white 
noise  receiver  was  achieved.   Low  noise  preamplifiers  used  to 
improve  receiver  sensitivity  are  also  a  source  of  colored 
noise.   The  effect  such  preamplifiers  have  on  practical  re- 
ceivers was  discussed  in  Chapter  V.   The  conclusion  of  such 
analysis  indicated  again  that  receivers  designed  to  operate 
in  a  colored  noise  environment  do  not  perform  significantly 
better  than  white  noise  receivers,  operating  in  a  similar 
environment.   Improvement  is  achieved  only  due  to  the  fact 
that  the  low  noise  preamplifier  isolates  the  front-end  of  the 
receiver  from  the  "noisy"  correlator. 

Although  these  results  were  obtained  using  a  first  order 
Butterworth  filter  amplifier  excited  by  white  noise  as  a 
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model  for  the  colored  noise  source,  it  is  reasonable  to  expect 
(based  on  the  results  of  Ref.  16)  that  using  a  more  compli- 
cated model  for  the  colored  noise  generation  will  not  change 
significantly  the  results  and  the  conclusions  obtained. 

The  ECM  vulnerability  of  the  colored  noise  receiver 
versus  that  of  the  white  noise  receiver  was  analyzed  in 
Chapter  VI.   Both  receivers  were  found  to  be  equally  sensi- 
tive to  a  deterministic  waveform  jammer.   However  when  the 
jammer  used  an  optimum  waveform  which  was  related  to  the 
signals  used  in  the  correlation  operation  in  the  receivers,  a 
significant  deterioration  in  the  performance  of  the  colored 
noise  receiver  was  observed.   It  must  be  noted  however  that 
the  jammer  waveform  designer  has  in  practice  almost  no  chance 
to  accurately  determine  what  this  optimum  waveform  should  be 
and  how  to  appropriately  use  it. 

This  thesis  has  demonstrated  the  relativ.e  robustness 
of  the  white  noise  receiver.   In  most  practical  cases,  even 
when  colored  noise  interferences  are  present,  the  white  noise 
receiver  performs  nearly  as  well  as  the  colored  noise  receiver 
designed  for  specific  interference  models.   Only  in  very 
special  cases  in  which  optimum  signal  waveforms  were  used  to 
transmit  the  binary  information  did  the  colored  noise  receiver 
perform  better  than  a  white  noise  receiver. 

Table  7.1  presents  a  quantitative  summary  of  the  perfor- 
mances of  the  colored  noise  receivers  and  white  noise  receivers 
Table  7.2  presents  the  effect  of  jammers  on  the  performances 
of  the  white  noise  and  the  colored  noise  receiver. 
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Table  7.1 


The  Receiver 

P 
e 

SNR  =  10,  JSR  =  1,  E 

=  1 

White  noise  receiver  rectangular 
pulses 

0.129 

Colored  noise  receiver  rectangular 
pulses 

0.125 

Colored  noise  receiver  sinusoidal 
pulses  b  =  1 ( *) 

0.07 

Colored  noise  receiver  sinusoidal 
pulses  b  =  4 

0.01 

Colored  noise  receiver  sinusoidal 
pulses  b  =  8 

0.00045 

White  noise  receiver  sinusoidal 
pulses  b  =  1 

0.28 

White  noise  receiver  sinusoidal 
pulses  b  =  4 

0.03 

White  noise  receiver  sinusoidal 
pulses  b  =  6 

0.002 

(*)  E  is  defined  by  Eq.  (4.44)  and  is  the  ratio  of  the  signal 
frequency  to  interference  noise  bandwidth. 


Table  7.2 


The  Receiver 

P 
e 

P 
e 

Jammer  waveform  equals 

Optimal  jammer 

signal  waveform 

SNR  =  10,  JSR  =  1 

SNR  =  10,  JSR  =  1,  E  =  0.1 

E  =  0.1 

JSR*  =  -5  db 

JSR*  =  -5  db 

White  noise  receiver 

rectangular  pulses 

0.05 

0.05 

Colored  noise  re- 

ceiver rectangular 

pulses 

0.05 

0,4 

*  JSR*  is  related  to  the  jammer  and  defined  as  E  ./E    ,    namely 

n~i   s 
the  ratio  of  interference  energy  to  signal      J 

energy. 
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APPENDIX  A 

DETAILED  SOLUTION  OF  A  FREDHOLM  II  EQUATION  FOR  COLORED 
NOISE  WITH  RATIONAL  SPECTRA  AND  RECTANGULAR  PULSES 

The  Fredholm  II  equation  to  be  solved  is 


N,  T 

-yhd(t)  +   /   Kc(t-u)hd(u)du   =   yd(t)   0£t<T     (A.l) 


The  noise  is  assumed  to  be  a  sample  function  from  a  W.S.S. 
process  whose  P.S.D.  is  given  by 


2 

<b  (s)   =   ~ j      7~  (A. 2) 

C         -s^  +3       D(sZ) 


The  autocorrelation  function  corresponding  to  this  P.S.D. 
is  given  by 


Kc(T)   =   cte  P|T|  (A. 3) 


The  signal  y^(t)  is  defined  by 


2A     0  <  t  <  T 


yd(t)   =  (A. 4) 

'  0      t  <  0,  t  >  T 


Equation  (A. 2)  can  be  written  as  follows 


D(s2)  [d>  (s)]   =   N(s2)   =   2a3  (A. 5) 

rc 
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2 

Multiplying  c|>  (s)  by  D(s  )  corresponds  to  operating  on  K  (t) 
*•*  c 

2 

with  D(p  )  in  the  time  domain  where  p  is  the  derivative 

operator. 

Eq.  (A. 5)  can  be  written  as 


D(p2) [Kc(t) ]   =   2a36(t)  (A. 6) 


2 

Operating  with  D(p  )  on  Eq.  (A.l)  yields 


2  N,  T     , 

d(p  )  -2"hd(t)  +   /  °tP  )  [Kc(t-cr)hd(a)]da  =  D(p2)[yd(t)]   (A. 7) 


Substituting  Eq.  (A. 6)  in  Eq.  (A.l)    and  performing  the 

2 
operation  specified  by  D(p  ),  yields  the  differential  equation 


N  ..        N  "2 

"TLhd(t)  +  ~T  6  hd(t)  +  2a^hd(t)   =  "  ^d(t)  +  B  Yd(t)    (A*8) 


Substituting  y,  (t)  as  defined  by  Eq.  (A. 4)  yields 


N 

-j   h{t)  +  6*2hit)   =   82A     0  <  t  <  T  (A. 9) 


where 


N 
2       1   2 
*    =  ~y   3   +  2a[ 


Eq.  (A. 9)  is  a  differential  equation.   Its  particular  solution 
is  given  by  a  constant  C,  where 
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3.2 

3' 


C   =   (^)z  2A  (A. 10) 


Its  homogeneous  solution  is 


K1eYt  +  K2e.  yt  (A. 11) 


where 


Y   =   /2/N,  3*  (A. 12) 


The  complete  solution  is  of  the  form 


h^t)   =   C  +  CK1eYt  +  CK2e  Yt  (A. 13) 


where  the  constants  K,  and  K„  are  obtained  by  plugging  this 
solution  into  the  Fredholm  II  equation. 

This  process  is  very  long  and  tedious  and  involves  a 
great  deal  of  algebraic  manipulations.   At  the  end  of  this 
process  two  linear  equations  are  obtained  which  define  K, 
and  K_ ,    namely 

"Kl  +   K2 


Y  +3    Y  -3     3 
K  pyT   k  p-yt 

Kl£      K2£         1 

Y"  3     y  +3      e 

The  simultaneous  solution  of  Eq.  (A. 14)  is 
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K. 


j(Y2  -e2)  [(y  -3)e    6T  +    (Y  +6)] 
"(y  +B)2eyT   -    (Y  -6)2e"yT 


K, 


12         2 
j(Y      -SZ)  [(Y 


)     +     (Y   +6)£YT] 


(Y  +B)2eyT 


(Y  -6)2e   yT 


where 


and 


Y      =      /2/N,    B*      = 


—  (-j.  3      +  2ctB) 
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(A. 15) 


c     =      (^r)2    2A     = 


2AB 


N 

-±B     +2aB 


(A. 16) 


=      B    /l  +NQ/N1 


(A. 17) 


APPENDIX  B 

DETAILED  SOLUTION  OF  A  FREDHOLM  II  EQUATION  FOR  COLORED 
NOISE  WITH  RATIONAL  SPECTRA  AND  SINUSOIDAL  PULSE  INPUT 

The  Fredholm  II  equation  to  be  solved  is 


N,  T 

-y  hd(t)  +   /   Kc(t-u)hd(u)du   =   yd(t)   0<t<T     (B.l) 


The  colored  noise  P.S.D.  and  autocorrelation  function  are 
the  same  as  defined  in  Appendix  A,  Eqs .  (A. 2),  (A. 3).   The 
signal  is  however  different  and  is  defined  as 

|  2A  sinbt     0  <  t  <_   T 
YA(t)  (B.2) 

0  >  t,  t  >  T 


I 


The  procedure  outlined  in  Appendix  A,  Eqs.  (A. 5) -(A. 8)  is 
applicable  here.   Using  Eq.  (A. 8)  and  using  y,(t)  as  derived 
by  Eq.  (B.2)  yields 


N 

-^h,(t)    +    3*2h.(t)       =      2A(b2  +32)sin  bt  (B.3) 

2d  d 


The  homogeneous  solution  is  not  affected  by  the  sinusoidal 
driving  function  in  this  equation.   Thus,  the  homogeneous 
solution  is  given  by  Eq.  (A. 11).   The  particular  solution  is 
given  by 
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h    (t)       =      C    sin  bt  (B.4) 


where  c    is    given  by 


C      =      N   2Alt>2  +&2]  (B.5) 

12  2 


The    complete    solution    is    therefore 


h(t)       =      C    sin   bt    +    CK^1    +    CK2e    Yt      0    <    t    <_  T       (B.6) 


d 


where  the  constants  K.  and  K»  are  obtained  by  plugging  this 
solution  into  the  Fredholm  II  equation.   This  process  is  very 
long  and  tedious  and  involves  a  great  deal  of  algebraic 
manipulations.   At  the  end  of  this  process  two  linear  equations 
are  obtained  which  define  K,  and  K? ,  namely 

"Kl  K2  -b 


y+e     y"e         e2+b2 


1  T1  ^9  ^  1 

eY      -   -=— ■ — =-     =        0         0  [B  sinbT  +b  cosbT] 


(B.7) 


Y-B         Y+  3      e2+b2 


The  solution  for  K,  and  K~  is 
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K. 


(|)  (*-!)£    2yT   +    (sin  bT  +|  cosbT)  (I+l)e"YT 


1-1 


(B.8) 


[1    +     (f)     ]  [-    H 


;-2YT 


+  1 


^+1 


i-i' 


K. 


-  |(1  +1)    +    (sin  bT  +  (|)cos   bT)  (J-1)e   YT 


2      1+1         1-1 
3  S 


(B.9) 


-2yT, 


129 


APPENDIX  C 
THE  FREDHOLM  II  EQUATION  FOR  BANDPASS  SIGNALS 

The  Fredholm  II  equation  for  bandpass  signals  and  spectra 
is  given  by 


N,  ~  T  ~ 

~Y   hd(t)  +   /   Kc(t-u)hd(u)du   =   yd(t)   0<t<T     (C.l) 


where 


y^,(t)   =  y,(t)  cos  co.t  (C.2) 


K  (t)   =   K  (t)  cos  wnt  (C.3) 

^  \*>  \J 


We  assume  that  the  solution  of  Eq.  (C.l)  is  of  the  form 


hd(t)   =  hd(t)  cos  ooQt  (C.4) 


and  check  the  conditions  under  which  this  assumption  is  valid. 
Substituting  Eqs.  (C.2),  (C.3),  and  (C.4)  into  Eq.  (C.l)  yields 


N,  T 

~2-hd(t)cos  wQt  +  j      K  (t-u)cos  to  (t-u)h,(u)cos  ooQudu 


=   yd(t)  cos  ajQt  (C.5) 
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Using  trigonometric  identities  yields  the  following  equation 


Nx  T 

-yhit)  cos  u)Qt  +  J      cos  co0tK  (t-u)hd(u)du 


T 
+   f   cos  oirt(t-2u)K  (t-u)h,  (u)du 
0         °        c       d 


=   Yd(t)  cos  o)Qt  (C.6) 


We  now  denote 


T 

/   cos  ojrt(t-2u)K  (t-u)h,(u)du     a(t)  (C.7) 

0         ° 


and  check  the  conditions  under  which  a(t)  is  negligible.   We 
may  write 


h*(t)   =   hd(t)  [u(t)  -  u(t-T)]  (C.8) 


Substituting  h^(t)  in  Eq.  (C.7)  enables  changing  the  limits 
of  integration 


K  (t-u)cos  oon  (t-2u)h*(u)du   =   a(t)  (C.9) 

c  u       a 


Taking  the  Fourier  transform  of  Eq.  (C.9)  yields 


00      oo 


F[a(t)]   =     /    /   K  (t-u)cos  co  (t-2u)h*(u)e'ja3tdudt     (CIO) 


c 

—00    — oo 
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with 


t-u      =      a  (C.ll) 


substituting   in   Eq.     (CIO)    yields 


F{a(t)}      =  /   h*(u)[      /    K    (a)cosa)Ja-u)e    3^°da]e   3amdu 

'       a  J       c  0 

—  oo  — oo 

(C.12) 

The  expression  inside  the  parenthesis  is  the  Fourier  transform 
of  K  (a)  cos  oon(a-u)  and  is  given  by  the  convolution  of  the 
Fourier  transform  of  K  (a)  and  cos  u>n(a-u).   Performing  the 
convolution  yields 


1  _:|aJ0U    1  JW0U 

7T   <j)  (oo  -coA)e       +  ~-  $  (a)  +u)n)e  (C.13) 

2  c      0  2   c      0 


where  <J>  (a))  is  the  P.S.D.  or  equivalently  the  Fourier  Transform 

of  K  (t)  . 
c 

Substituting   Eq.     (C.13)     into    Eq .     (C.12)    yields 


F[a(t)]       =      ±        j      h*(u)  (cj)c  (  go-co  0)e         U 


—  00 


]0)nU 

+    (j>c(a)  +oo0)£       Q    )£_:WUdu  (C.14) 


Evaluating    the    integral   yields 


^■H*(oo  +600)cj)c(a)  -co0)     +    |   4>c(co  +coQ)H*(co  -L0Q)  (C.15) 

132 


Eq.  (C.15)  is  a  cross  product  of  two  low  pass  functions  at 
high  frequencies.   If 


0       T 


then  Eq.  (C.15)  represents  a  negligible  small  value.   Under 
these  conditions,  a(t)  is  a  small  number. 

The  bandpass  Eq.  (C.6)  can  be  written  as 

N,  T 

-T-h,(t)    cosoint   +   cos   oo~t      f      K    (t-u)h,(u)du      =      y,(t)cos   cont 
2.    d  0  O^c  d  Jd  0 

(C.16) 

The  solution  to  Eq.  (C.16)  is  the  solution  to  the  lowpass 
equation  multiplied  by  cos  wnt. 

The  conclusion  is  that  if  the  communication  channel 
center  frequency  is  much  bigger  than  the  channel  bandwidth, 
then  the  solution  to  Eq.  (C.l)  can  be  written  as  a  solution 
to  a  baseband  equation  multiplied  by  cos  cont  where  con  is 
the  channel  center  frequency. 
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APPENDIX  D 

BLOCK  DIAGRAM  OF  A  SIGNAL  GENERATOR  FOR  h,(t) 
a 


A  suggested  block  diagram  for  generating  h,  (t)  in  the 
correlator  unit  is  shown  in  Fig.  D.l. 

The  desired  waveform  is  sampled  at  a  high  enough  rate. 
The  samples  are  digitized  and  stored  in  the  PROM.   When  the 
sync  pulse  is  received,  the  proper  addresses  of  the  PROM  are 
read  sequentially  and  the  output  is  converted  to  an  analog 
signal  h,  (t) .   The  PROM  can  store  several  waveforms  for 
various  types  of  transmission  signals  or  interferences. 
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APPENDIX  E 
DETAILED  SOLUTION  OF  THE  FREDHOLM  II  EQUATION  GIVEN  BY  EQ.  (5.9) 

The  Fredholm  II  equation  to  be  solved  is 

N  T 

-|-hd(t)  +   /   Kc(t-u)hd(u)du   =   2AG(l-e"3t)   =   yd(t)   (E.l) 

0  £  t  <_   T 
where 


Kc(T)   =   ae-SM   =  ^BG2  «-<»M  (E.2) 


The  solution  follows  the  procedure  shown  in  Appendix  A 
except  that  y,  is  no  longer  a  constant  and  therefore  modifica- 
tion to  the  particular  solution  will  be  necessary.   The 
procedure  of  Appendix  A  is  applicable  here  up  to  Eq.  (A. 8), 
so  that  our  starting  point  is 


N0  "       N0  2  "2 

-f  h^t)  +  ^-s  ^(t)  +  2a3h^t)   =   -yd(t)  +  6  yd(t)        (E.3) 


Substituting  yd(t)  as  defined  by  Eq.  (E.l)  yields 


N   ..         ?  2 

-j-  hit)  +  8*  hit)   =   2AS  G  (E.4) 


where 
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2      N0   2  N0   2    N0   2  2      N0  2      2 

3*z   =   _£.  g^  +  2a3   =   -|  r  +  -j  TG^   =   -j-B  (1  +  G^)    (E.5) 


The  homogeneous  solution  to  Eq.  (E.4)  is 


CK1eyt  +  CK2e  yt  (E.6) 


where  y    is  given  by 


Y   =   /2/Nq  3*   =   3  A  +G2   =   3mx 


(E.7) 


The  particular  solution  is  given  by 

2 

C   =   (JL)  2AG   =   4AG  «  (E.8) 

P  NQ(1  +GZ) 

The  complete  solution  therefore  is 


hd(t)   =  C  +  CK1eyt   +  CK2e'Yt 


The  next  step  is  to  plug  this  solution  into  Eq.  (E.l) .   This 
substitution  leads  to  two  equations  that  must  be  solved  for 
K,  ,  and  K~  ,  that  is, 

~K1  +   K2        1 


Y  +3    Y  -6       3 

(E.9) 

K  pyt   „  -yt 

Kl£      K2£         1 

Y-3~Y+3     3 
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The  simultaneous  solution  to  these  equations  is 


£-[(Y  +6)  -  (y  -3)e  YT]  (Y2  -32) 


Kl   =  YT,  ,ft,2 J -,2  -YT  (E'10) 

e '  (y  +  8)   -  (y  -  8)  e  ' 


i(Y2  -62)  t(Y  -6)  -  (Y  +6)e  TT] 

K2  =  " 1  ~2   YT J -.2  -YT  (E-11^ 

(Y  +$)  e  '   -  (y  -8)  £  ' 


Substituting  Eq.  (E.7)  into  Eqs.  (E.10)  and  (E.ll),  and 
introducing  the  notation 

E   =   8T 

yields 


-Em    m, -1   -2Em. 

(mn+l) [e    X  -  -ixv  £    X] 

l  m,  +1 

K    = (E.12) 

m,  +1    m. -1   -2Em. 

m, -1    m. +1 

-Em.    m  +1 

(m.-l)  [e    L   -   -^r] 
1  m.  -1 

K2   = (E.13) 

m. +1    m. -1   -2Em. 
m, -1    m. +1 
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